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1.1 Introduction 

Chiral symmetries play an important role in the spectrum and phenomenology of both 
the standard model and various theories for physics beyond the standard model. In 
many cases chiral symmetry is associated with nonperturbative physics which can only 
be quantitatively explored in full on a lattice. It is therefore important to implement 
chiral symmetry on the lattice, which turns out to be less than straightforward. In 
these lectures I discuss what chiral symmetry is, why it is important, how it is broken, 
and ways to implement it on the lattice. There have been many hundreds of papers on 
the subject and this is not an exhaustive review; the limited choice of topics I cover 
reflects on the scope of my own understanding and not the value of the omitted work. 

1.2 Spinor representations of the Lorentz group 

To understand chiral symmetry one must understand Lorentz symmetry first. Since 
we will be discussing fermions in various dimensions of spacctime, consider the gener- 
alization of the usual Lorentz group to d dimensions. The Lorentz group is defined by 
the real matrices A which preserve the form of the d-dimensional metric 

A T V A = r,, r ? = diag(l,-l,...-l). (1.1) 

With this definition, the inner product between two 4-vectors, v^rj^w^ — v T r]w, is 
preserved under the Lorentz transformations v — > Av and w — > Aw. This defines 
the group SO(d — 1, 1), which — like SO(d) — has d(d — l)/2 linearly independent 
generators, which may be written as M^ v = —M v ^, where the indices /i, v — 0, . . . , (d— 
1) and 

A = e i6 ^ M ^ , (1.2) 

with 6u V — —Qvij, being d(d— l)/2 real parameters. Note that fi, v label the d(d — l)/2 
generators, while in a representation R each M is a da x da matrix, where da is the 



dimension of R. By expanding eqn. (1.1) to order 9 one sees that the generators M 
must satisfy 

(M^frj + rjMf" = 0. (1.3) 

From this equation it is straightforward to write down a basis for the M^ v in the d- 
dimensional defining representation and determine the commutation relations for the 
algebra, 
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[M af3 , M^] = i (rj^M 06 - -q^M 05 - n p *M cr ' + r] aS ' M^'') . (1.4) 
A Dirac spinor representation can be constructed as 

M a ? = S a/3 = l - [ 7 Q , V*] (1.5) 

where the gamma matrices satisfy the Clifford algebra: 

{ 7 Q , = Irf? (1.6) 

Solutions to the Clifford algebra are easy to find by making use of direct products 
of Pauli matrices. In a direct product space we can write a matrix ds M — a® A where 
a and A are matrices of dimension d a and d& respectively, acting in different spaces; 
the matrix M then has dimension (d a x <1a)- Matrix multiplication is defined as (a (g> 
A) (6® B) = (a&) (8) (-AS). It is usually much easier to construct a representation when 
you need one rather than to look one up and try to keep the conventions straight! One 
finds that solutions for the 7 matrices in <i-dimcnsions obey the following properties: 

1 . For both d = 2k and d = 2k + 1 , the 7-matrices are 2 k dimensional; 

2. For even spacetime dimension d = 2k (such as our own with k = 2) one can define 
a generalization of 75 to be 

2fe-l 

r = i k - 1 Yl 7 M (i-7) 

with the properties 

{r, 7 ^} = o, r = r t = r- 1 , Tr(r 7 Qi • •• 7 Q2fc ) = 2 k r 1 - k e ai - a2k , (1.8) 

where e i2...2fe-i = +1 = -e 012 " 2fe - 1 . 

3. In d — 2k + 1 dimensions one needs one more 7-matrix than in d = 2k, and one 
can take it to be 7 2fe = iT. 

Sometimes it is useful to work in a specific basis for the 7-matrices; a particulary useful 
choice is a "chiral basis", defined to be one where T is diagonal. For example, for d = 2 
and d = 4 (Minkowski spacetime) one can choose 

d = 2 : 7 = en , 7 1 = -ic7 2 , r = cr 3 (1.9) 

rf = 4: 7° = -(Ti®l, 7 4 = ier 2 ® cr, , T = cr 3 (g) 1 . (1.10) 

1.2.1 7-matrices in Euclidian spacetime 

In going to Euclidian spacetime with metric i}^ = 5^ u , one takes 

d™^id*, df^df (1.11) 



and defines 
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J°m = ^e, lii = i!%, (1.12) 

so that 

(7^=7^, {7^,7l} = 2^ (1.13) 

and Tp M — > with 

$ £ = -lf E (1.14) 

and the Euclidian Dirac operator is + m). The matrix L( 2fe ) in 2k dimensions is 
taken to equal 7^ in (2fc + 1) dimensions: 

li 2fc) - 7f = r£ fe) , Tr(r B 7* ■ ■ - 7 £") = -2 fc i*e ai - ' a * (1.15) 
where e i 2 ...2fc-i = +1 = + e 012 - 2fe - 1 . 

1.3 Chirality in even dimensions 

My lectures on chiral fermions follow from the properties of the matrix T in even di- 
mensions. The existence of T means that Dirac spinors are reducible representations of 
the Lorentz group, which in turn means we can have symmetries ("chiral symmetries") 
which transform different parts of Dirac spinors in different ways. To see this, define 
the projection operators 

P ± = <^ , (1.16) 

which have the properties 

P+ + P- = l, P^=P±, P+P-=0. (1.17) 
Since in odd spatial dimensions {L,7 M } = for all /i, it immediately follows that V 



commutes with the Lorentz generators Yi^ v in eqn. (1.5|: [r, E M "] = 0. Therefore we 



can write W = Yf£ + TT where 

^ = P ± E^P ± , E^E^ = E^E^f = . (1.18) 

Thus E Ml/ is reducible: spinors ij)± which are eigenstates of L with eigenvalue ±1 
respectively transform independently under Lorentz transformations. 

The word "chiral" comes from the Greek word for hand, x eL P- The projection 
operators P + and P_ are often called Pr and Pl respectively; what does handedness 
have to do with the matrix L? Consider the Lagrangian for a free massive Dirac fcrmion 
in 1 + 1 dimensions and use the definition of L = 7°7 1 : 
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= & + i (d t + d x ) + #tt (d t - d x ) (1.19) 

where 

r*± = ±*± . (1.20) 
We sec then that the solutions to the Dirac equation for m — are 

9±(x,t) = f±(iTi) 

so that 9 + corresponds to right-moving solutions, and \E'_ corresponds to left-moving 
solutions. This is possible since the massless particles move at the speed of light and 
the direction of motion is invariant under proper Lorentz transformations in (1 + 1) 
dimensions. 

In the chiral basis eqn. (1.9 1, the positive energy plane wave solutions to the Dirac 
equation are 

= e -^(*— ) Q , y_ = e -is(t+x) M (1 21) 

It is natural to call P+ and P- as Pr and Pl respectively. 

Exercise 1.1 You should perform the same exercise in 3 + 1 dimensions and find that so- 
lutions ty± to the massless Dirac equation satisfying Y^± = ±\l/± must also satisfy \p\ — E 
and (2p • S/E)*Sf± — ±$±, where Si = ^eoijh^ are the generators of rotations. Thus i Sf± 
correspond to states with positive or negative helicity respectively, and are called right- and 
left-handed particles. 



1.4 Chiral symmetry and fermion mass in four dimensions 

Consider the Lagrangian for a single flavor of Dirac fermion in 3+1 dimensions, coupled 
to a background gauge field 

C = 9(0 - m)* = (V l 09l + 9r09r) - m (# L tf fl + * R V L ) (1.22) 

where I have defined 

* L = P_*, * L = *[ 7 = , ¥ H = P + tf , * fl = *P_. (1.23) 

For now I am assuming that 9 l.r arc in the same complex representation of the gauge 
group, where is the gauge covariant derivative appropriate for that representation. 
It is important to note the the property {75, 7^} = ensured that the kinetic terms in 



eqn. (1.231 do not couple left-handed and right-handed fcrmions; on the other hand, 
the mass terms dcP] The above Lagrangian has an exact U(l) symmetry, associated 
with fermion number, * —* e la \D '. Under this symmetry, left-handed and right-handed 
components of \& rotate with the same phase; this is often called a "vector symmetry" . 



1 I will use the familiar 75 in 3 + 1 dimensions instead of T when there is no risk of ambiguity. 
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In the case where m = 0, it apparently has an additional symmetry where the left- and 
right-handed components rotate with the opposite phase, — > e 1 " 75 ^; this is called 
an "axial symmetry", U(1)a- 

Symmetries are associated with Noether currents, and symmetry violation appears 
as a nonzero divergence for the current. Recall the Noether formula for a field (f> and 
infinitesimal transformation d> —> <t> + cSd: 



J» = - a/a ^ 5(j) , d„J fi = -8£. (1.24) 



dC 

In the Dirac theory, the vector symmetry corresponds to = , and the axial 
symmetry transformation is d^f — 275 'J, so that the Noether formula yields the vector 
and axial currents: 

[7(1) : J" = , d^J" =0 (1.25) 

U(l) A : J>X = *y7s* , 9 M J£ =2im*7 5 tf . (1.26) 

Some comments arc in order: 



• Eqn. ( 1.26 1 is not the whole story! We will soon talk about additional contributions 
to the divergence of the axial current from the regulator, called anomalies, which 
do not decouple as the regulator is removed. 

• The fact that the fermion mass breaks chiral symmetry means that fermion masses 
get multiplicatively renormalized, which means that fcrmions can naturally be 
light (unlike scalars in most theories); more on this later. 

• The variation of a general fermion bilinear ^X^ under chiral symmetry is 

S^X^f = W {75, X}* . (1.27) 

This will vanish if X can be written as the product of an odd number of gamma 
matrices. In any even dimension the chirally invariant bilinears include currents, 
with X — 7^ or X = 7 M r, while the bilinears which transform nontrivially under 
the chiral symmetry include mass terms, X = 1, T. Thus gauge interactions can be 
invariant under chiral symmetry, while fermion masses are always chiral symmetry 
violating. In (3 + 1) dimensions, anomalous electromagnetic moment operators 
corresponding to X — u^ v , 07^75 are also chiral symmetry violating. 

• A more general expression for the classical divergence of the axial current for a 
bilinear action J ^fD^> in any even dimension is 

d^ = i^{T,D}^ . (1.28) 

On the lattice one encounters versions of the fermion operator D which violate 
chiral symmetry even for a massless fermion. 

The Lagrangian for Nf flavors of massive Dirac fermions in odd d, coupled to some 
background gauge field may be written as 

l = (nupn + - (nM ab * b R + ^m^i) . (1.29) 

The index on ^> denotes flavor, with a,b = 1, . . . Nf, and M a f, is a general complex mass 
matrix (no distinction between upper and lower flavor indices). Again assuming the 
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fermions to be in a complex representation of the gauge group, this theory is invariant 
under independent chiral transformations if the mass matrix vanishes: 

*%->U ab * b R , * a L ^V ab V b L , rtU = V*V = l. (1.30) 

where U and V are independent U(Nf) matrices. Since U(Nf) = SU(Nf) x U(l), it 
is convenient to write 

U = e^ a+ ^R, V = e i( - a ~^L, tfR = tfL = l, \R\ = \L\ = 1 , (1.31) 

so that the symmetry group is SU(N f ) L xSU(N f ) R xU(l)xU(l) A with L G SU(N f ) L , 

ReSU(N f ) R . 

If we turn on the mass matrix, the chiral symmetry is explicitly broken, since the 
mass matrix couples left- and right-handed fermions to each other. If M a b = mb a b 
then the "diagonal" or "vector" symmetry SU(Nf) x U(l) remains unbroken, where 



SU(Nf) C SU(Nf)z x SU(Nf)n corresponding to the transformation eqn. (1.30), 
eqn. ( 1.31| with L = R. If M a b is diagonal but with unequal eigenvalues, the sym- 



metry may be broken down as far as U(l) Nf , corresponding to independent phase 
transformations of the individual flavors. With additional flavor-dependent interac- 
tions, these symmetries may be broken as well. 

1.5 Weyl fermions 

1.5.1 Lorentz group as SU(2) x SU(2) and Weyl fermions 

We have seen that Dirac fermions in even dimensions form a reducible representation 
of the Lorentz group. Dirac notation is convenient when both LH and RH parts of 
the Dirac spinor transform as the same complex representation under a gauge group, 
and when there is a conserved fermion number. This sounds restrictive, but applies 
to QED and QCD. For other applications — such as chiral gauge theories (where 
LH and RH fermions carry different gauge charges, as under SU(2) x U(l)), or when 
fermion number is violated (as is the case for neutrinos with a Majorana mass), or 
when fermions transform as a real representation of gauge group — then it is much 
more convenient to use irreducible fermion representations, called Weyl fermions. 

The six generators of the Lorentz group may be chosen to be the three Hermitian 
generators of rotations Ji, and the three anti- Hermitian generators of boosts Ki, so 
that an arbitrary Lorentz transformation takes the form 

A = eWiJi+^K*) . (i.32) 



In terms of the M^ v generators in \1.2 



Ji = \^uM^ , Ki = M 0i . (1.33) 
These generators have the commutation relations 

[Ji,Jj\ = itijkJk , [Ji,Kj] = ie ijk K k , [K it Kj] = -ie i]k J k . (1-34) 



It is convenient to define different linear combinations of generators 
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M = , Bl = J -^i , (1.35) 

satisfying an algebra that looks like SU(2) x SU(2), except for the fact that the 



generators eqn. (1.351 are not Hermitian and therefore the group is noncompact: 



[A h A,} = ie ijk A k , [Bi,Bj] = ie ljk B k , [A t , Bj] = . (1.36) 

Thus Lorentz representations may be labelled with two SU{2) spins, ]a,b correspond- 
ing to the two SU(2)s: {ja, 3b), transforming as 

A(0, uj) = D jA (8- iuo) x D jB (6 + uj) (1.37) 

where the Z?- 7 is the usual SU{2) rotation in the spin j representation; boosts appear 
as imaginary parts to the rotation angle; the D 3A and D JB matrices act in different 
spaces and therefore commute. For example, under a general Lorentz transformation, 
a LH Weyl fermion tp and a RH Weyl fermion \ transform as tp — > Lip, \ ~ ¥ R-X-> 
where 

L = e l{S - 1 ^- 3 ' 2 , R = e^^/ 2 . (1.38) 

Evidently the two types of fermions transform the same way under rotations, but 
differently under boosts. 

The dimension of the {ja, ]b) representation is (2j J 4 + l)(2js + l)- In this notation, 
the smaller irreducible Lorentz representations are labelled as: 

(0,0) : scalar 

(|, 0) , (0, \) : LH and RH Weyl fermions 

(|, i) : four- vector 

(1,0) , (0, 1) : self-dual and anti-self-dual antisymmetric tensors 

A Dirac fermion is the reducible representation (|,0) © (0, \) consisting of a LH and 
a RH Weyl fermion. 

Parity interchanges the two SU (2)s, transforming a (ji, j'2) representation into 



(j2, ji)- Similarly, charge conjugation effectively flips the sign of Ki in eqn. (1.371 due 
to the factor of i implying that if a field transforms as (ji, j'2), then (jy transforms as 
(j2; ii J] Therefore a theory of Nl flavors of LH Weyl fermions ipi, and Nr flavors of 
RH Weyl fermions \a ma Y be recast as a theory of (Nl + Nr) LH fermions by defining 
Xa = <^>a- The fermion content of the theory can be described entirely in therms of LH 
Weyl fermions then, {tpi, cu a }', this often simplifies the discussion of parity violating 
theories, such as the Standard Model or Grand Unified Theories. Note that if the RH 
Xa transformed under a gauge group as representation R, the conjugate fermions ui a 
transform under the conjugate representation R. 



2 For this reason, the combined symmetry CP does not alter the particle content of a chiral theory, 
so that CP violation muast arise from complex coupling constants. 
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For example, QCD written in terms of Dirac fermions has the Lagrangian: 

C= *n{0-m n )$ n , (1-39) 



i—u^d.s. 



where £)„ is the SU(3) C covariant derivative, and the ^ n fields (both LH and RH 
components) transform as a 3 of SU(3) C . However, we could just as well write the 
theory in terms of the LH quark fields ip n and the LH anti-quark fields Xn- Using the 



7-matrix basis in eqn. ( 1.10 1, we write the Dirac spinor in terms of two-component 
LH spinors tp and \ as 

Note that ip transforms as a 3 of SU(3) C , while \ transforms as a 3. Then the kinetic 
operator becomes (up to a total derivative) 

= tpUD^ip + x^iD^x , = {1, -3} , (1-41) 

and the mass terms become 

*_r*l = X^ = ip<T2X 

^L^B = ^<J 2 X f = X^2^ , (1-42) 

where I used the fact that fermion fields anti-commute. Thus a Dirac mass in terms 
of Weyl fermions is just 

= m(ipa 2 x + h.c) , (1-43) 

and preserves a fermion number symmetry where if) has charge +1 and x nas charge 
— 1. On the other hand, one can also write down a Lorentz invariant mass term of the 
form 

m(-!/w 2 V> + h.c.) (1-44) 

which violates fermion number by two units; this is a Majorana mass, which is clumsy 
to write in Dirac notation. Experimentalists are trying to find out which form neutrino 
masses have — Dirac, or Majorana? If the latter, lepton number is violated by two 
units and could show up in neutrinoless double beta decay, where a nucleus decays by 
emitting two electrons and no anti-neutrinos. 

The Standard Model is a relevant example of a chiral gauge theory. Written in 
terms of LH Weyl fermions, the quantum numbers of a single family under SU(3) x 
SU(2) x 17(1) are: 

Q = (3,2) + i i=(l,2)_| 
C/ C = (3,l)_| £ c = (l,l) +1 

7J c = (3,l) + i. (1.45) 

Evidently this is a complex representation and chiral. If neutrino masses are found 
to be Dirac in nature (i.e. lepton number preserving) then a partner for the neutrino 
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must be added to the theory, the "right handed neutrino" , which can be described 
by a LH Weyl fermion which is neutral under all Standard Model gauge interactions, 
JV=(l,l)o. 

If unfamiliar with two-component notation, you can find all the details in Ap- 



pendix A of Wess and Bagger's classic book on supersymmetry (Wess and Bagger, 



1992 1; the notation used here differs slightly as I use the metric and 7-matrix conven- 



tions of Itzykson and Zuber (Itzykson and Zuber, 19801, and write out the 172 matrices 
explicitly. 



Exercise 1.2 Consider a theory of Nf flavors of Dirac fermions in a real or pseudo-real 
representation of some gauge group. (Real representations combine symmetrically to form an 
invariant, such as a triplet of SU(2); pseudo-real representations combine anti-symmetrically, 
such as a doublet of SU(2)). Show that if the fermions are massless the action exhibits 
a U(2Nf) = 17(1) x SU(2Nf) flavor symmetry at the classical level (the U(l) subgroup 
being anomalous in the quantum theory). If the fermions condense as in QCD, what is the 
symmetry breaking pattern? How do the resultant Goldstone bosons transform under the 
unbroken subgroup of SU(2Nf)l 



Exercise 1.3 To see how the (|, |) representati on beh aves like a four-vector, consider the 
2x2 matrix P — P A ,cr M , where <j m is given in eqn. ( 1.41 1. Show that the transformation P — » 
LPL+ (with det L = 1) preserves the Lorentz invariant inner product P^P M = (P 2 
Show that with L given by eqn. (1.38 1, P M transforms properly like a four- vector. 



PP) 



Exercise 1.4 Is it possible to write down an anomalous electric or magnetic moment oper- 
ator in a theory of a single charge-neutral Weyl fermion? 



1.6 Chiral symmetry and mass renormalization 

Some operators in a Lagrangian suffer from additive renormalizations, such as the unit 
operator (cosmological constant) and scalar mass terms, such as the Higgs mass in the 
Standard Model, \H\ 2 . Therefore, the mass scales associated with such operators will 
naturally be somewhere near the UV cutoff of the theory, unless the bare couplings 
of the theory are fine-tuned to cancel radiative corrections. Such fine tuning problems 
have obsessed particle theorists since the work of Wilson and 't Hooft on renormaliza- 
tion and naturalness in the 1970s. However, such intemperate behavior will not occur 
for operators which violate a symmetry respected by the rest of the theory: if the bare 
couplings for such operators were set to zero, the symmetry would ensure they could 
not be generated radiatively in perturbation theory. Fermion mass operators generally 
fall into this benign category. 

Consider the following toy model: QED with a charge-neutral complex scalar field 
coupled to the electron: 

£f (i# - m)tf + \d(b\ 2 - n 2 \4>? - 3I0I 4 + V (*b0*£ + #i>**h) • (1-46) 

Note that in the limit m — > this Lagrangian respects a chiral symmetry W — > e 1 " 75 ^, 
4> — > e~ 2la tfi. The symmetry ensures that if m — 0, a mass term for the fermion would 
not be generated radiatively in perturbation theory. With m ^ 0, this means that 
any renormalization of m must be proportional to m itself (i.e. m is "multiplicatively 
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Fig. 1.1 One-loop renormalization of the electron mass in QED due to photon exchange. A 
mass operator flips chirality, while gauge interactions do not. A contribution to the electron 
mass requires an odd number of chirality flips, and so there has to be at least one insertion 
of the electron mass in the diagram: the electron mass is multiplicatively renormalized. A 
scalar interaction flips chirality when the scalar is emitted, and flips it back when the scalar 
is absorbed, so replacing the photon with a scalar in the above graph again requires a fermion 
mass insertion to contribute to mass renormalization. 




Fig. 1.2 One-loop additive renormalization of the scalar mass due to a quadratically diver- 
gent fermion loop. 

renormalized"). This is evident if one traces chirality through the Fcynman diagrams; 
see Fig. |1.1| Multiplicative renormalization implies that the fermion mass can at most 
depend logarithmically on the cutoff (by dimensional analysis): 5m ~ (a/47r)m In m/A. 
In contrast, the scalar mass operator \4>\ 2 does not violate any symmetry and 



therefore suffers from additive renormalizations, such as through the graph in Fig. 1.2 
By dimensional analysis, the scalar mass operator can have a coefficient that scales 
quadratically with the cutoff: 5/i 2 ~ (y 2 /16n 2 )A 2 . This is called an additive renormal- 
ization, since S/i 2 is not proportional to ll 2 . It is only possible in general to have a 
scalar in the spectrum of this theory with mass much lighter than yA/Ait if the bare 
couplings arc finely tuned to cause large radiative corrections to cancel. When referring 
to the Higgs mass in the Standard Model, this is called the hierarchy problem. 

If chiral symmetry is broken by operators other than the mass term, then the 
fermion mass will no longer in general be multiplicatively renormalized, and fine tun- 
ing may be necessary. This is particularly true if chiral symmetry is broken by "ir- 
relevant" operators. Consider adding to QED a dimension five operator of the form 
W — D |U -D^ 1 E r , where r is a dimensionless coupling. This operator breaks chirality 
and therefore one can substitute it for the mass operator in Fig. |1.1| an estimate of the 
diagram then gives an additive renormalization of the fermion mass, 5m ~ (a/47r)r , A. 
Thus unless r is extremely small (e.g. r < m/A) the chiral symmetry breaking effects 
of this operator will be important and fine tuning will be necessary to ensure a light 
fermion in the spectrum. This example is relevant to Wilson's method for putting 
fcrmions on the lattice, which does not respect chiral symmetry and entails adding 
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to the action a lattice version of W, where a ~ 1/A is the lattice spacing and r ~ 1. 
Therefore Wilson fermions acquire an 0(1/ a) correction to their mass which needs to 
be canceled by a bare contribution in order to describe a world with light fermions. 

1.7 Chiral symmetry in QCD 

1.7.1 Chiral symmetry breaking and Goldstone bosons 

So far the discussion of chiral symmetry in terms of the effect on fermion masses has 
been appropriate for a weakly coupled theory. As was presented in M. Golterman's 
lectures, the low energy spectrum of QCD is described by a chiral Lagrangian, encod- 
ing the interactions of the meson octet which are the approximate Nambu- Goldstone 
bosons of spontaneously broken SU(3)l X SU(3)r chiral symmetry]^] The Goldstone 
bosons would be massless in the limit of exact chiral symmetry, and so tuning away 
the leading finite lattice space correction for Wilson fermions can be accomplished by 
tuning the bare quark mass to eliminate the 1 /a dependence of the square of the pion 
mass. 

In contrast to QCD, N — 1 super Yang Mills theory has a single Weyl fermion (the 
gaugino) transforming as an adjoint under the gauge group; the theory has a U(1)a 
symmetry at the classical level — phase rotations of the gaugino — but it is broken 
by anomalies to a discrete symmetry. This discrete symmetry is then spontaneously 
broken by a gluino condensate, but without any continuous symmetries, no Goldstone 
bosons are produced. What should the spectrum of this theory look like? Presumably 
a bunch of massive boson and fermion glucball-like states. They will form degenerate 
supersymmetric multiplets when the gluino mass is tuned to zero, but there is no 
particle that becomes massless in the chiral limit in this case, and therefore tuning the 
bare mass is difficult. 

After tuning away the 0(1/ a) mass correction, there remain for non-chiral lattice 
fermions the dimension-5 chiral symmetry violating operators in the Symanzik action 
which require 0(a) tuning, as discussed by Golterman. In contrast, chiral fermions 
receive finite lattice corrections only at 0(a 2 ), simply because one cannot write down 
a dimcnsion-5 chiral symmetry preserving operator in QCD. 

1.7.2 Operator mixing 

One encounters additional factors of 1/a when computing weak processes. One of the 
most curious feature of the strong interactions is the AI = 1/2 rule, which is the 
observation that As = 1 transitions in nature are greatly enhanced when they change 
isospin by AI = 1/2, in comparison to AI = 3/2. For example, one requires for the 
amplitudes for kaon decay K — > tttt: 

3 With six favors of quarks in QCD, one might ask why an SU(3) chiral symmetry instead of 5(7(2) 
or SU (6). The point is that the chiral symmetry is broken by quark masses, and whether the breaking 
is large or small depends on the ratio m q / Aqqd, where here Aqcd is some strong interaction scale 
in the 100s of MeV. The u and d quarks are much lighter than Aq^jj, the strange quark is borderline, 
and the c, b, t quarks are much heavier. Therefore SU(2) X 5(7(2) is a very good symmetry of QCD; 
5(7(3) X 5(7(3) is a pretty good symmetry of QCD, but assuming chiral symmetry for the heavier 
quarks is not justified. Radiative corrections in the baryon sector go as \Jra q j Kqq-q and so there 
even SU (3) X SU (3) does not appear to be very reliable. 
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A(AI = 3/2)~ M - [ 

To compute this in the standard model, one starts with four-quark operators generated 
by W-exchange, which can be written as the linear combination of two operators 

C AS =i = -VudV:^ [C+{n,M w )0+ + C-(n,M„)0-] , 

± = [{sd) L {uu) L ± (su) L (ud) L ] -[u^c] , (1.48) 

where {qq')L = (97^ Plq')- If onc ignores the charm quark contribution, the 0~ trans- 
forms as an 8 under SU(3)f, while + transforms as a 27; therefore 0~ is pure 
I = 1/2, while + is a mix of I = 3/2 and 7 = 1/2. The full ± operators are in 
SU{A)f multiplets; while SU(4)j is not a good symmetry of the spectrum, it is only 
broken by quark masses which do not effect the log divergences of the theory. Thus 
the running of the operators respect SU(A)f down to /i = m c , and there is no mixing 
between ± . At the weak scale fi = Myy, one finds |C + /C~| = 1 + 0(a s (Mw), show- 
ing that there is no AI — 1/2 enhancement intrinsic to the weak interactions. One 
then scales these operators down to fi ~ 2 GeV in order to match onto the lattice the- 
ory; using the renormalization group to sum up leading a s In ^/Mw corrections gives 
an enhancement |C + /C~| ~ 2 — which is in the right direction, but not enough to 



explain eqn. (1.47), which should then cither be coming from QCD at long distances, 
or else new physics! This is a great problem for the lattice to resolve. 

A wonderful feature about using dimensional regularization and MS in the con- 
tinuum is that an operator will never mix with another operator of lower dimension. 
This is because there is no UV mass scale in the scheme which can make up for the 
miss-match in operator dimension. This is not true on the lattice, where powers of the 
inverse lattice spacing X/a can appear. In particular, the the dimension-6 four fermion 
operators ± could in principle mix with dimension-3 two fermion operators. The 
only AS = 1 dimension-3 operator that could arise is s^d, which is also AI — \\\ 
If the quarks were massless, the lattice theory would possess an exact discrete "SCP 1 ^ 
symmetry under which on interchanges s <-> d and performs a CP transformation to 
change LH quarks into LH anti-quarks; the operators 0± arc even under SCP while 
s"/§d is odd, to the operator that could mix on the lattice is 

O p = (m s ~ m d ) sj 5 d . (1.49) 

In a theory where the quark masses are the only source of chiral symmetry breaking, 
then O p = d^A^f, the divergence of the AS = 1 axial current. Therefore on-shell 
matrix elements of this operator vanish, since the derivative gives (px — P2tt) — 0, 
i.e. no momentum is being injected by the weak interaction. We can ignore O p then 
when the K — > ttt: amplitude is measured with chiral lattice fermions with on-shell 
momenta. 



4 The operator sd is removed by rediagonalizing the quark mass matrix and does not give rise to 
Kitty. 
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For a lattice theory without chiral symmetry, O p — d^Aff + 0(a) and so has a 



nonvanishing 0(a) matrix clement. In this case operators 0± from eqn. (1.481 in the 
continuum match onto the lattice operators 



± {o) + %O p 



0(a) . (1.50) 



In general then one would need to determine the coefficient C p to O(a) in order to 
determine the AI = | amplitude for K — » irir to leading order in an a-expansion, 
which is not really feasible. Other weak matrix elements such as Bk and e'/e similarly 
benefit from the use of lattice fermions with good chiral symmetry. 

1.8 Fermion determinants in Euclidian space 

Lattice computations employ Monte Carlo integration, which requires a positive inte- 
grand that can be interpreted as a probability distribution. While not sufficient for a 
lattice action to yield a positive measure, it is certainly necessary for the continuum 
theory one is approximating to have this property. Luckily, the fermion determinant 
for vector- like gauge theories (such as QCD), det(j0 + m), has this property in Euclid- 
ian space. Since 1$ = —ij) and {T.lfi} = 0, it follows that there exist eigenstates ip n 
of Jp such that 

Ifiipn = iXnipn , IpTtp n = ~i\ n Tip n , X n real. (1-51) 

For nonzero A, ip m and Tip n are all mutually orthogonal and we see that the eigenvalue 
spectrum contains ±iA n pairs. On the other hand, if A„ = then ip n can be an 
eigenstate of T as well, and Tip n is not an independent mode. Therefore 

det($ + m)= l[l(\l + m 2 )x JJ m (1.52) 

A„>0 A rl =0 

which is real and for positive m is positive for all gauge fields. 

What about a chiral gauge theory? The fermion Lagrangian for a LH Wcyl fermion 



in Euclidian space looks like ipD^ip with — D^a^ and (in the chiral basis eqn. (1.10), 
continued to Euclidian space) = {l,i(?}. Note that Dl has no nice hermiticity 
properties, which means its determinant will be complex, its right eigenvectors and 
left eigenvectors will be different, and its eigenvectors will not be mutually orthogo- 
nal. Furthermore, Dl is an operator which maps vectors from the space C of LH Weyl 
fermions to the space 1Z of RH Wcyl fermions. In Euclidian space, these spaces are 
unrelated and transform independently under the SU(2) x SU(2) Lorentz transfor- 
mations. Suppose we have an orthonormal basis \n,TZ) for the RH Hilbert space and 
\n,C) for the LH Hilbert space; we can expand our fermion integration variables as 

4> = y^c ra ln.,£) , ijj = y]cn{n,ll\ (1.53) 



so that 
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Fig. 1.3 The eigenvalue How of the Dirac operator as a function of gauge fields, and two 
unsatisfactory ways to define the Weyl fermion determinant det Dl as a square root of detlfi. 
The expression | det lj> | corresponds to the picture on the left, where det Dl is defined as the 
product of positive eigenvalues of Ij); this definition is nonanalytic at A t . The picture on the 
right corresponds to the product of half the eigenvalues, following those which were positive 
at some reference gauge held Ao . This definition is analytic, but not necessarily local. Both 
definitions are gauge invariant, which is incorrect for an anomalous fermion representation. 

/[#][#] e-^ 73 ^ =det(m,K\D L \n,£) . (1.54) 

However, the answer we get will depend on the basis we choose. For example, we could 
have chosen a different orthonormal basis for the C space \n',L) = U n > n \n,C) which 
differed from the first by a unitary transformation U ; the resultant determinant would 
differ by a factor detli, which is a phase. If this phase were a number, it would not 
be an issue — but it can in general be a functional of the background gauge field, so 
that different choices of phase for det Dl lead to completely different theories. 

We do know that if Dr is the fermion operator for RH Weyl fermions in the same 
gauge representation as Dl, then detD,R = detD* L and that det Dr det Dl = Act Ifi. 
Therefore the the norm of |det-Di| can be defined as 

det D L = ^\detIp\(? w[A] ( 1 . 55) 

where the phase W[A] is a functional of the gauge fields. What do we know about 
W[A]1 

1. Since det0 is gauge invariant, W\A] should be gauge invariant unless the fermion 
representation has a gauge anomaly, in which case it should correctly reproduce 
that anomaly; 

2. It should be analytic in the gauge fields, so that the computation of gauge field 
correlators (or the gauge current) are well defined. 

3. It should be a local functional of the gauge fields. 

In Fig. 1 1 . 3 1 1 show two possible ways to define det Dl, neither of which satisfy the 
above criteria. The naive choice of just setting W[A] — not only fails to reproduce 
the anomaly (if the fermion representation is anomalous) but is also nonanalytic and 
nonlocal. It corresponds to taking the product of all the positive eigenvalues A n of ifi 
(up to an uninteresting overall constant phase). This definition is seen to be nonan- 
alytic where eigenvalues cross zero. Another definition might be to take the product 
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of positive eigenvalues at some reference gauge field A$, following those eigenvalues 
as they cross zero; this definition is analytic, but presumably not local, and is always 
gauge invariant. 

There has been quite a few papers on how to proceed in defining this phase W[A] 
in the context of domain wall fermions, including a rather complicated explicit con- 
struction for U(l) chiral gauge theories on the lattice ( Luscher, 1999| [Luscher, 20 00a); 
however, even if a satisfactory definition of W[A] is devised, it could be impossible to 
simulate using Monte Carlo algorithms due to the complexity of the fermion determi- 
nant. 



1.9 Parity and fermion mass in odd dimensions 

In these lectures I will be discussing fermions in (2k + 1) dimensions with a spatially 
varying mass term which vanishes in some 2fc-dimensional region; in such cases we 
find chiral modes of a 2fc-dimensional effective theory bound to this mass defect. Such 
an example could arise dynamically when fermions have a Yukawa coupling to a real 
scalar cf> which spontaneously breaks a discrete symmetry, where the surface with tj> = 
forms a domain wall between two different phases; for this reason such fermions are 
called domain wall fermions, even though we will be putting the spatially dependent 
mass in by hand and not through spontaneous symmetry breaking. 

To study domain wall fermions it is useful to say a few words about fermions in 
odd dimensions where there is no analogue of T and therefore there is no such thing 
as chiral symmetry. Nevertheless, fermion masses still break a symmetry: parity. In a 
theory with parity symmetry one has extended the Lorentz group to include improper 
rotations: spatial rotations R for which the determinant of R is negative. Parity can 
be defined as a transformation where an odd number of the spatial coordinates flip 
sign. In even dimensions parity can be the transformation x — ► — x and 

#(x,£) -> 7°*(-x,t) (parity, d even) . (1.56) 

Note that under this transformation and are exchanged and that a Dirac mass 
term is parity invariant. 

However, in odd dimensions the transformation x — > — x is just a rotation; instead 
we can define parity as the transformation which just flips the sign of one coordinate 
x , and 

#(x,i) -> 7 1 *(x,t) , x= (~x\x 2 ,...,x 2k ) (1.57) 

Remarkably, a Dirac mass term flips sign under parity in this case; and since there 
is no chiral symmetry in odd d to rotate the phase of the mass matrix, the sign of 
the quark mass is physical, and a parity invariant theory of massive quarks must have 
them come in pairs with masses ±M, with parity interchanging the two. 



1.10 Fermion masses and regulators 

We have seen that theories of fermions in any dimension can possess symmetries which 
forbid masses - chiral symmetry in even dimensions and parity in odd dimensions. 
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This property obviously can have a dramatic impact on the spectrum of a theory. 
Supersymmetry ingeniously puts fermions and bosons in the same supermultiplct, 
which allows scalars to also enjoy the benefits of chiral symmetry, which is one reason 
theorists have been so interested in having supersymmetry explain why the Higgs boson 
of the standard model manages to be so much lighter than the Planck scale. However, 
precisely because mass terms violate these symmetries, it is difficult to maintain them 
in a regulated theory. After all, one regulates a theory by introducing a high mass 
scale in order to eliminate UV degrees of freedom in the theory, and this mass scale 
will typically violate chiral symmetry in even dimensions, or parity in odd. This gives 
rise to "anomalous" violation of the classical fcrmion symmetries, my next topic. 



2 

Anomalies 



2.1 The U(1)a anomaly in 1+1 dimensions 

One of the fascinating features of chiral symmetry is that sometimes it is not a sym- 
metry of the quantum field theory even when it is a symmetry of the Lagrangian. In 
particular, Noether's theorem can be modified in a theory with an infinite number of 
degrees of freedom; the modification is called "an anomaly" . Anomalies turn out to 
be very relevant both for phenomenology, and for the implementation of lattice field 
theory. The reason anomalies affect chiral symmetries is that regularization requires a 
cut-off on the infinite number of modes above some mass scale, while chiral symmetry 
is incompatible with fermion massed 

Anomalies can be seen in many different ways. I think the most physical is to 
look at what happens to the ground state of a theory with a single flavor of massless 
Dirac fermion in (1 + 1) dimensions in the presence of an electric field. Suppose one 
adiabatically turns on a constant positive electric field E(t), then later turns it off; the 
equation of motion for the fermion is [2] ^ = e E(t) and the total change in momentum 
is 



Ap=e J E{t)dt . (2.1) 



Thus the momenta of both left- and right-moving modes increase; if one starts in the 
ground state of the theory with filled Dirac sea, after the electric field has turned off, 
both the right-moving and left-moving sea has shifted to the right as in Fig. |2.1| The 
the final state differs from the original by the creation of particle- antiparticle pairs: 
right moving particles and left moving antiparticles. Thus while there is a fermion 
current in the final state, fermion number has not changed. This is what one would 
expect from conservation of the U(l) current: 

V" = , (2.2) 
However, recall that right-moving and left-moving particles have positive and negative 



chirality respectively; therefore the final state in Fig. 2.1 has net axial charge, even 
though the initial state did not. This is peculiar, since the coupling of the electromag- 
netic field in the Lagrangian does not violate chirality. We can quantify the effect: if 

1 Dimensional regularization is not a loophole, since chiral symmetry cannot be analytically con- 
tinued away from odd space dimensions. 

2 While in much of these lectures I will normalize gauge fields so that D M = c9 M + iA^, in this 
section I need to put the gauge coupling back in. If you want to return to the nicer normalization, 
set the gauge coupling to unity, and put a 1/g 2 factor in front of the gauge action. 
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Fig. 2.1 On the left: the ground state for a theory of a single massless Dirac fermion in 
(1 + 1) dimensions; on the right: the theory after application of an adiabatic electric field 
with all states shifted to the right by Ap, given in eqn. ( |2.1[ ). Filled states are indicated by 
the heavier blue lines. 

we place the system in a box of size L with periodic boundary conditions, momenta 
are quantized as p n = 2irn/L. The change in axial charge is then 



AQ A = 2 



Ap 
2%/L 



d 2 xE(t) = 



2tt 



d 2 x e 



//// 



pttu 



(2.3) 



where I expressed the electric field in terms of the field strength F, where F 01 = 
— F 10 = E. This can be converted into the local equation using AQa = J cPxd^J 1 ^, a 
modification of eqn. (1.26): 



<9 M J*X = 2im^T^> 



2tt 



(2.4) 



where in the above equation I have included the classical violation due to a mass term 
as well. The second term is the axial anomaly in 1 + 1 dimensions; it would vanish for 
a nonabelian gauge field, due to the trace over the gauge generator. 

So how did an electric field end up violating chiral charge? Note that this analysis 
relied on the Dirac sea being infinitely deep. If there had been a finite number of 
negative energy states, then they would have shifted to higher momentum, but there 
would have been no change in the axial charge. With an infinite number of degrees of 
freedom, though, one can have a "Hilbert Hotel" : the infinite hotel which can always 
accommodate another visitor, even when full, by moving each guest to the next room 
and thereby opening up a room for the newcomer. This should tell you that it will not 
be straightforward to represent chiral symmetry on the lattice: a lattice field theory 
approximates quantum field theory with a finite number of degrees of freedom — the 
lattice is a big hotel, but quite conventional. In such a hotel there can be no anomaly. 

We can derive the anomaly in other ways, such as by computing the anomaly dia- 



gram Fig. 2.2 or by following Fujikawa (Fujikawa, 1979 Fujikawa, 1980) and carefully 



accounting for the Jacobian from the measure of the path integral when performing 
a chiral transformation. It is particularly instructive for our later discussion of lattice 
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Fig. 2.2 The anomaly diagram in 1+1 dimensions, with one Pauli-Villars loop and an in- 
sertion of 2iM$r$ at the X. 

fermions to compute the anomaly in perturbation theory using Pauli-Villars regulators 
of mass M. We replace our axial current by a regulated current 

J A,re g = #VT# + 57T$ , (2.5) 

where <& is our Pauli-Villars field; it follows then that 

d n J A,reg = 2im*r* + 2iM5 r$ . (2.6) 

We are interested in matrix elements of rcg in a background gauge field between 
states without any Pauli-Villars particles, and so we need to evaluate (2iM<I>r<I>) in 
a background gauge field and take the limit M — > oo to see if 9 M picks up any 

anomalous contributions that do not decouple as we remove the cutoff. 

To compute (2iM$r$) we need to consider all Feynman diagrams with a Pauli- 
Villars loop, and insertion of the <f>r<f> operator, and any number of external U(l) 
gauge fields. By gauge invariance, a graph with n external photon lines will contribute 
n powers of the field strength tensor F^ v . For power counting, it is convenient that 
we normalize the gauge field so that the covariant derivative is — (c^ + iA^); then 
the gauge field has mass dimension 1, and has dimension 2. In (1 + 1) dimen- 
sions (2iM $r<i>) has dimension 2, and so simple dimensional analysis implies that the 
graph with n photon lines must make a contribution proportional to (F P v ) n / M 2 ^ 1 ^ 1 ' . 
Therefore only the graph in Fig. |2.2| with one photon insertion can make a contribution 
that survives the M — > oo limit (the graph with zero photons vanishes). Calculation 
of this diagram yields the same result for the divergence of the regulated axial current 



as we found in eqn. (2.4 1 



Exercise 2.1 Compute the diagram in Fig. |2.2| using the conventional normalization of the 
gauge field D M = (<9 M + ieA^) and verify that 2iM ($r$) = ^e M „F'" / when M -> oo. 



Note that in this description of the anomaly we (i) effectively rendered the number 
of degrees of freedom finite by introducing the regulator; (ii) the regulator explicitly 
broke the chiral symmetry; (iii) as the regulator was removed, the symmetry break- 
ing effects of the regulator never decoupled, indicating that the anomaly arises when 
the two vertices in Fig. 2.2 sit at the same spacetime point. While we used a Pauli- 
Villars regulator here, the use of a lattice regulator will have qualitatively similar 
features, with the inverse lattice spacing playing the role of the Pauli-Villars mass, 
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and we can turn these observations around: A lattice theory will not correctly repro- 
duce anomalous symmetry currents in the continuum limit, unless that symmetry is 
broken explicitly by the lattice regulator. This means we would be foolish to expect 
to construct a lattice theory with exact chiral symmetry. But can the lattice break 
chiral symmetry just enough to explain the anomaly, without losing the important 
consequences of chiral symmetry at long distances (such as protecting fermion masses 
from renormalization)? 



2.2 Anomalies in 3+1 dimensions 

2.2.1 The U(l) A anomaly 

An analogous violation of the U(1)a current occurs in 3 + 1 dimensions as well |^] 
One might guess that the analogue of e^F^ = 2E in the anomalous divergence 



eqn. (2.4 1 would be the quantity e Alz , po .F A " y F p<T = 8E-B, which has the right dimensions 
and properties under parity and time reversal. So we should consider the behavior a 
massless Dirac fermion in (3 + 1) in parallel constant E and B fields. First turn on a 
B field pointing in the z direction: this gives rise to Landau levels, with energy levels 
E n characterized by non-negative integers n as well as spin in the z direction S z and 
momentum p z , where 

E l =P 2 z + (2" + l)eB - 2eBS z . (2.7) 

The number density of modes per unit transverse area is defined to be g ni which 
can be derived by computing the zero-point energy in Landau modes and requiring 
that it yields the free fermion result as B — ► 0. We have g n — ► p±dp±/(2ir) with 
[(2n + l)eB - 2eBS z ] -> p 2 x , implying that 

g n = eB/2n . (2.8) 

Th dispersion relation looks like that of an infinite number of one-dimensional fermions 
of mass m n ,±, where 

m 2 n ± = (2n+l)eB-2eBS z , S z = ±\ . (2.9) 

The state with n — and S z = +5 is distinguished by having m n ,+ — 0; it behaves 
like a massless one-dimensional Dirac fermion (with transverse density of states go) 
moving along the i axis with dispersion relation E — \p z \ . If we now turn on an electric 
field also pointing along the z direction we know what to expect from our analysis in 
1 + 1 dimensions: we find an anomalous divergence of the axial current equal to 



goeE/n = e'EB/2n' - (^^J W^^ P " ■ ( 2 -!0) 
If we include an ordinary mass term in the 3 + 1 dimensional theory, then we get 



3 Part of the content of this section comes directly from John Preskill's class notes on the strong 
interactions, available at his web page: http://www.theory.caltech.edu/'~preskill/notes.html. 



Anomalies in 3+1 dimensions 21 




Fig. 2.3 The U(1)a anomaly diagram in 3+1 dimensions, with one Pauli-Villars loop and 
an insertion of 2iAf<E>r$. 



= 2im^T^ 



167r 2 ; 



(2.11) 



One can derive this result by computing (MQiTQ) for a Pauli-Villars regulator as 
in the 1 + 1 dimensional example; now the relevant graph is the triangle diagram of 
Fig.lO 



If the external fields are nonabelian, the analogue of eqn. (2.11 1 is 

-,2 



JL \ f pP a Tr T 71 

^g 7r 2 / t VP<? r a r b ±L -La-Lb 



(2.12) 



If the fermions transform in the defining representation of SU (N) , it is conventional 
to normalize the coupling g so that Tr T a Tb = ^S a i,. This is still called an "Abelian 
anomaly", since generates a 17(1) symmetry. 

2.2.2 Anomalies in Euclidian spacetime 



Continuing to Euclidian spacetime by means of eqns. ( |l.ll[ )- pTT5 ) changes the anomaly 
equations simply by eliminating the factor of i from in front of the fermion mass: 



2d : d^ = 2m^r* + -e^f" 

Z7T 



( 9 2 



Ad : d^ A = 2m«T* + ( ^ ) e^FfFg* Tr T a T b 



V16tt 5 



2.2.3 The index theorem in four dimensions 



(2.13) 
(2.14) 



For nonabelian gauge theories the quantity on the far right of eqn. (2.14) is a topolog- 
ical charge density with 



•L f d i Tvf VPV-ppcr 

64w 2 J ■ CEt ^P" r a r a 



(2.15) 



being the winding number associated with ir$(G), the homotopy group of maps of £3 
(spacetime infinity) into the gauge group G. 



Consider then continuing the anomaly equation eqn. (2.12) to Euclidian space and 



integrating over spacetime its vacuum expectation value in a background gauge field 
(assuming the fermions to be in the TV-dimensional representation of SU(N) so that 
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Tr T a Tb = \5 a b)- The integral of 9 M (J^) vanishes because it is a pure divergence, so 
we get 



d 4 x E m(*r*) = -v . (2.16) 
The matrix element above on the right equals 

J [d«f}[dV] er SE (m#r#) I J [d^][d^] e- SE . (2.17) 

where Se = ^{IPe + Tn)^. We can expand * and ^ in terms of eigenstates of the 
anti-hermitian operator J^ E , where 

lJ> E ^n = i^nlpn , J d i X E ^ m ^n = 8 mn , (2-18) 

with 



Then 

d 4 XE m 



m)= (]T I d 4 x E miPiTiP n l[(i\ k + m)) /\{{i\ k + m) 

= mJ2 f d 4 x E ipnTip n /(i\i + m) . (2.20) 

n 

Recall that {T,$} = 0; thus 

$ipn = iK^n implies $(Tip n ) = -i\ n (Tip n ) . (2.21) 

Thus for X n 7^ 0, the eigenstates ip n and (r^) must be orthogonal to each other (they 
are both eigenstates of if) with different eigenvalues), and so iplTip n vanishes for A„ =^ 
and does not contribute to the sum in eqn. ( |2.20 1. In contrast, modes with A„ = can 



simultaneously be eigenstates of lj> and of T; let n + , n_ be the number of RH and LH 
zeromodes respectively The last integral in then just equals (n + — n_) = {tir — Ul), 



and combining with eqn. (2.16) we arrive at the index equation 

n- - n+ = v , (2.22) 

which states that the difference in the number of LH and RH zeromode solutions to 
the Euclidian Dirac equation in a background gauge field equals the winding number 
of the gauge field. With Nf flavors, the index equation is trivially modified to read 

n_ - n + = NfV . (2.23) 

This link between eigenvalues of the Dirac operator and the topological winding num- 
ber of the gauge field provides a precise definition for the topological winding number 
of a gauge field on the lattice (where there is no topology) — provided we have a 
definition of a lattice Dirac operator which exhibits exact zeromodes. We will see that 
the overlap operator is such an operator. 
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Fig. 2.4 Anomalous three-point function of three currents. 
2.2.4 More general anomalies 

Even more generally, one can consider the 3-point correlation function of three arbi- 
trary currents as in Fig. |2.4| 



mk).f b {p)J2{q)) 



(2.24) 



and show that the divergence with respect to any of the indices is proportional to a 
particular group theory factor 



M J£(A0 -TO Jf(g)) « Tr Qa{Qb, Qc} 



R—L 



(2.25) 



where the Qs are the generators associated with the three currents in the fcrmion 
representation, the symmetrized trace being computed as the difference between the 
contributions from RH and LH fcrmions in the theory. The anomaly A for the fcrmion 
representation is defined by the group theory factor 



Tr {Q a {Qb,Qc}) 



R-L 



Ad a 



be •) 



(2.26) 



with d a bc being the totally symmetric invariant tensor of the symmetry group. For 
a simple group G (implying G is not U(l) and has no factor subgoups), d a b c is only 
nonzero for G = SU(N) with N > 3; even in the case of SU(N), d a bc will vanish 
for real irreducible representations (for which Q a = —Q*), or for judiciously chosen 
reducible complex representations, such as 5 © 10 in SU{5). For a semi-simple group 
Gi x Gi (where G\ and G^ are themselves simple) there are no mixed anomalies 
since the generators are all traceless, implying that if Q € G\ and Q € Gi then 
Tr (Q a {Qb, Qc}) oc Tr Q a = 0. When considering groups with U(l) factors there 
can be nonzero mixed anomalies of the form U(l)G 2 and U(l) 3 where G is simple; 
the U(l) 3 anomalies can involve different U(l) groups. With a little group theory it 
is not difficult to compute the contribution to the anomaly of any particular group 
representation. 

If a current with an anomalous divergence is gauged, then the theory does not make 
sense. That is because the divergenceless of the current is required for the unphysical 
modes in the gauge field to decouple; if they do not decouple, their propagator 
has a piece that goes as k^k u /k 2 which does not fall off at large momentum, and the 
theory is not renormalizablc. 
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When global U(l) currents have anomalous divergences, that is interesting. We 
have seen that the U(1)a current is anomalous, which explains the rf mass; the di- 
vergence of the axial isospin current explains the decay ir° — > 77; the anomalous 
divergence of the baron number current in background SU (2) in the Standard Model 
predicts baryon violation in the early universe and the possibility of weak-scale baryo- 
gencsis. 

Exercise 2.2 Verify that all t he ga uge currents are anomaly-free in the standard model 
with the representation in eqn. ( 1.45 1. The only possible G 3 anomalies are for G = SU(3) or 



G — U(l); for the SU(S) 3 anomaly use the fact that a LH Weyl fermion contributes +1 to 
A if it transforms as a 3 of SU(3), and contributes — 1 to A if it is a 3. There are two mixed 
anomalies to check as well: U(1)SU(2) 2 and U(1)SU(3) 2 . 

This apparently miraculous cancellation is suggestive that each family of fermions may be 
unified into a spinor of SO (10), since the vanishing of anomalies which happens automatically 
in 50(10) is of course maintained when the symmetry is broken to a smaller subgroup, such 
as the Standard Model. 

Exercise 2.3 Show that the global B ( baryo n number) and L (lepton number) currents are 



anomalous in the Standard Model eqn. ( 1.45 1, but that B — L is not 



2.3 Strongly coupled chiral gauge theories 

Strongly coupled chiral gauge theories are particularly intriguing, since they can con- 
tain light composite fermions, which could possibly describe the quarks and leptons 
we see. A nice toy example of a strongly coupled chiral gauge theory is SU(5) with 
LH fermions 

= 5 , X = 10 ■ (2.27) 
It so happens that the ip an d the x contribute with opposite signs to the (SU(5)) 3 



anomaly A in eqn. (2.26), so this seems to be a well defined gauge theory. Further- 
more, the SU(5) gauge interactions are asymptotically free, meaning that interactions 
becomes strong at long distances. One might therefore expect the theory to confine as 
QCD does. However, unlike QCD, there are no gauge invariant fermion bilinear con- 
densates which could form, and which in QCD are responsible for baryon masses. That 
being the case, might there be any massless composite fermions in the spectrum of 
this theory? 't Hooft came up with a nice general argument involving global anomalies 
that suggests there will be. 

In principle there are two global U(l) chiral symmetries in this theory correspond- 
ing to independent phase rotations for ip and x; however both of these rotations have 
global x SU(5) 2 anomalies, similar to the global x SU(3) 2 of the U(1)a current in QCD. 
This anomaly can only break one linear combination of the two U(l) symmetries, and 
one can choose the orthogonal linear combination which is anomaly-free. With a little 
group theory you can show that the anomaly-free global U(l) symmetry corresponds 
to assigning charges 

^ = 5 3 , X = 10-i, (2.28) 

where the subscript gives the global U(l) charge. This theory has a nontrivial global 
U(l) 3 anomaly, A = 5 x (3) 3 + 10 x (-1) 3 = 125. 't Hooft's argument is that 
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this(global) 3 anomaly restricts — and helps predict — the low energy spectrum of the 
theory. Applied to the present model, his argument goes as follows: imagine weakly 
gauging this U(l) symmetry. This would be bad news as the theory stands, since a 
(gauge) anomaly leads to a sick theory, but one can add a LH "spectator fermion" 
uj = 1_5 which is a singlet under SU(5) but has charge —5 under this U(l) symmetry, 
canceling the U(l) 3 anomaly. This weak U(l) gauge interaction plus the SU(5)-singlet 
uj fermion should not interfere with the strong SU{b) dynamics. If that dynamics 
leads to confinement and no U(l) symmetry breaking, then the weak U(l) gauge the- 
ory must remain anomaly free at low energy, implying that there has to be one or more 
massless composite fermions to cancel the U(l) 3 anomaly of the uj. A good candidate 
massless composite LH fermion is (tp^X) which is an S'?7(5)-singlet (as required by 
confinement), and which has U(l) charge of (3 + 3 — 1) = 5, exactly canceling the 
J7(l) 3 anomaly of the uj. Now forget the thought experiment: do not gauge the U(X) 
and do not include the u spectator fermion. It should still be true that this SU(5) 
gauge theory produces a single massless composite fermion (tpipx) Q 

While it is hard to pin down the spectrum of general strongly coupled chiral gauge 
theories using 't Hooft's anomaly matching condition alone, a lot is known about 
strongly coupled supersymmetric chiral gauge theories, and they typically have a very 
interesting spectrum of massless composite fermions, which can be given small masses 
and approximate the quarks and leptons we see by tweaking the theory. See for example 



(Kaplan, Lepeintre and Schmaltz, 19971, which constructs a theory with three families 
of massless composite fermions, each with a different number of constituents. 

Chiral gauge theories would be very interesting to study on the lattice, but pose 
theoretical problems that have not been solved yet — and which if they were, might 
then be followed by challenging practical problems related to complex path integral 
measures and massless fermions. Perhaps one of you will crack this interesting problem. 

2.4 The non-decoupling of parity violation in odd dimensions 

We have seen that chiral symmetry does not exist in odd space dimensions, but that a 
discrete parity symmetry can forbid a fermion mass. One would then expect a regulator 
- such as Pauli-Villars fields — to break parity. Indeed they do: on integrating the 
Pauli-Villars field out of the theory, one is left with a Chern Simons term in the 
Lagrangian with coefficient M/\M\, which does not decouple as M — ► oo. In 2k + 1 
dimensions the Chern Simons form for an Abelian gauge field is proportional to 

c iai J a 2 a 3 * a 2 ka2k + l 

which violates parity; the Chern Simons form for nonabelian gauge fields is more 
complicated. 

For domain wall fermions we will be interested in a closely related but slightly 
different problem: the generation of a Chern Simons operator on integrating out a 

4 You may wonder about whether fermion condensates form which break the global U(l) symmetry. 
Perhaps, but it seems unlikely. The lowest dimension gauge invariant fermion condensates involve four 
fermion fields — such as (xXXV") or UxV'Xx 1 /')^} — which are all neutral under the U(l) symmetry. 
Furthermore, there are arguments that a Higgs phase would not be distinguishable from a confining 
phase for this theory. 
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Fig. 2.5 Integrating out a heavy fermion in three dimensions gives rise to the Chern Simons 
term in the effective action of eqn. (2.301. 



heavy fermion of mass m. In 1 + 1 dimensions with an Abelian gauge field one computes 
the graph in Fig. |2.5| which gives rise to the Lagrangian 



2 

87r m 



(2.30) 



What is interesting is that it implies a particle number current 



_ I dCcs _ e m 



e 8tt |m| 



a/3 



(2.31) 



which we will see is related to the anomaly eqn. (2.4 1 in 1 + 1 dimensions. 



Exercise 2.4 Verify the coefficient in eqn. (2.301 by computing the diagram Fig. 2.5 By 
isolating the part that is proportional to e M ^ a p before performing the integral, one can 
make the diagram very easy to compute. 



3 

Domain Wall Fermions 



3.1 Chirality, anomalies and fermion doubling 

You have heard of the Nielsen-Ninomiya theorem: it states that a fermion action in 
2k Euclidian spacetime dimensions 

fir /a J2k _ 

S= -^*- p D(pMp) (3.1) 

cannot have the operator D satisfy all four of the following conditions simultaneously: 

1. D(p) is a periodic, analytic function of p^; 

2. D(p) oc 7 M p p for a\p^\ < 1; 

3. -D(p) invcrtible everywhere except p^ = 0; 

4. {r,£(p)} = o. 

The first condition is required for locality of the Fourier transform of D(p) in 
coordinate space. The next two state that we want a single flavor of conventional Dirac 
fermion in the continuum limit. The last item is the statement of chiral symmetry. One 
can try keeping that and eliminating one or more of the other conditions; for example, 
the SLAC derivative took D(p) — j^p^ within the Brillouin zone (BZ), which violates 
the first condition — if taken to be periodic, it is discontinuous at the edge of the BZ. 
This causes problems — for example, the QED Ward identity states that the photon 
vertex T M is proportional to dD(p)/dp fJ- , which is infinite at the BZ boundary. Naive 
fermions satisfy all the conditions except (3): there D(p) vanishes at the 2 4 corners 
of the BZ, and so we have 2 4 flavors of Dirac fermions in the continuum. Staggered 
fermions are somewhat less redundant, producing four flavors in the continuum for 
each lattice field; Creutz fermions are the least redundant, giving rise to two copies 
for each lattice field. The discussion in any even spacetime dimension is analogous. 

This roadblock in developing a lattice theory with chirality is obviously impossible 
to get around when you consider anomalies. Remember that anomalies do occur in 
the continuum but that in a UV cutoff on the number of degrees of freedom, there 
are no anomalies, and the exact symmetries of the regulated action are the exact 
symmetries of the quantum theory. The only way a symmetry current can have a 
nonzero divergence is if either the original action or the UV regulator explicitly violate 
that symmetry. The implication for lattice fermions is that any symmetry that is exact 
on the lattice will be exact in the continuum limit, while any symmetry anomalous in 
the continuum limit must be broken explicitly on the lattice. 
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A simple example to analyze is the case of a "naive" lattice action for a single RH 
fermion, 

S= 2^ f+(n)7 M [* + (n+A)-*+(n-A)] 
= ^H 2isina ^*+(-P)7A l *+(p) , 

r#+ = (3.2) 

so that D(p) = sin ak^/a. This vanishes at every corner of the BZ; expanding 
about these points we write = + n^ir/a with n M 6 {0, 1} and ag M -C 1 and find 

D(p)~^(-ir-y^ . (3.3) 

These zeroes of D(p) (at g M = 0) correspond to the 2 d doublers in d-dimensional Eu- 
clidian spacetime, violating condition (3) in the Nielsen-Ninomiya theorem. However, 
D(jp) does satisfy condition (4) and the action S is invariant under the symmetry 

*+(p) - e iQ *+(p) = e iar *+(p) , (3.4) 

which looks like a chiral symmetry — yet for a continuum theory with 2 d RH fermions, 
a phase symmetry would be anomalous, which we know cannot result from a symmetric 
lattice theory! 

The resolution is that the continuum theory does not have 2 d RH fermions, but 
rather 2 d ~ 1 Dirac fermions, and the exact lattice U(l) symmetry corresponds to an 
exact fermion number symmetry in the continuum, which is not chiral and not anoma- 



lous. To show this, note that D(p) in eqn. (3.3 1 has funny signs near the corners of 
the BZ. We can convert back to our standard gamma matrix basis using the similar- 
ity transformation P(n) ltJi P{n)- 1 = (-1)"^; but then P(ri)rP(n)- 1 = (-l) CT (")r, 
where o~(n) — n M (since V is the product of all the 7 M ). Therefore 

r[P(n)*+(q)] = (-ir<")[P(n)v]/ + (q)] (3.5) 

and in the continuum we have 2 d_1 RH fermions and 2 d ~ 1 LH fermions, and the exact 
and apparently chiral symmetry of the lattice corresponds to an exact and anomaly- 
free fermion number symmetry in the continuum. The redundancy of staggered and 
Creutz fermions serve the same purpose, ensuring that all exact lattice symmetries 
become anomaly-free vector symmetries in the continuum limit. 



3.2 Domain wall fermions in the continuum 

3.2.1 Motivation 

What we would like is a realization of chiral symmetry on the lattice which (i) is not 
exact, so we can correctly recover anomalies, but which (ii) retains all the good features 
of chiral symmetry in the continuum, such as protection from additive renormalization 
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of fermion masses. There is a curious example in the continuum of such a system, 
which gave a clue on how to achieve this. The example has to do with fermions in odd 
dimension that interact with a domain wall. To be concrete, consider a system in three 
dimensions (coordinates (xo, X\, #2)), where the fermion has a mass which depends on 
X2 and switches sign at x-i = 0. For simplicity, I will take mix-i) = mtix-i) = mxij\x-2\. 
Curiously enough, we will show that a massless fermion mode exists bound to this 2- 
dimcnsional surface, and that it is chiral: there exist a RH mode, and not a LH one. 
Thus the low energy limit of this theory looks like a 2-dimensional theory of a Weyl 
fermion with a chiral symmetry, even though we started with a 3-dimensional theory 
in which there can be no chirality. 

Yet we know that the low energy effective theory is anomalous. Recall that for a 
massless Dirac fermion coupled to photons in two Euclidian dimensions, the vector 
current is conserved ( d^J^ = 0), while the axial current is not in general {p^J^ = 
(e 2 /wE)). Thus the fermion current for a RH Weyl current satisfies 



«=^W) = ^- (3-6) 

If we turn on an electric field pointing in the X\ direction, then the charge on the mass 
defect must increase with time. Yet in the full 3-dimensional theory, there is only one 
fermion current J p = vp^i]/ , and it is conserved. So even though only massive states 
live off the mass defect, we see that they must somehow know about the anomaly and 
allow chiral symmetry to be violated as the anomaly requires. 

You might be suspicious that there is some hidden fine tuning here to keep the 
chiral mode massless, but that cannot be: the low energy effective theory would need 
a LH mode as well in order for there to be a mass. Distortions of the domain wall 
mass function cannot change this result, unless the mass m(x2) becomes small enough 
somewhere to change the spectrum of the low energy effective theory. But even in that 
case there is an index theorem that requires there to be a massless Weyl fermion as 
the sign of m{x2) changes at an odd number of locations. 

This looks like a useful trick to apply to the lattice: an anomalous chiral symmetry 
emerging at low energy from a full theory with no fundamental chiral symmetry, and 
without fine tuning. In this lecture I show how the continuum theory works, and then 
how it can be transcribed to the lattice. In the next lecture I will discuss how the 
effective theory can be described directly using the overlap formulation, without any 
reference to the higher dimensional parent theory. 

3.2.2 The model 

Even though fermions in even and odd dimensions look quite different, one finds an 
interesting connection between them when considering the Dirac equation with a space- 
dependent mass term. One can think of a space-dependent mass as arising from a Higgs 
mechanism, for example, where there is a topological defect trapped in the classical 
Higgs field, such as a domain wall or a vortex. A domain wall can naturally arise when 
the Higgs field breaks a discrete symmetry; a vortex when the Higgs field breaks a 
U(l) symmetry (see John Preskill's lectures "Vortices and Monopoles" at the 1985 



Les Houches Summer School (Preskill, 1985). Domain wall defects are pertinent to 



putting chiral fermions on the lattice, so I will consider that example. 
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Consider a fermion in Euclidian spacetime with dimension d — 2k + 1 , where the 
coordinates are written as {xq, x%,... X2k-i, s} = s}, where fx = 0, . . . , 2k— 1 and s 
is what I call the coordinate x^k- The (2/c+l) 7 matrices are written as {70, . . . , 72k, T}. 
This fermion is assumed to have an s-dependent mass with the simple form 



m(s) = me(s) 



+m s > 
—to s < 



to > 



(3.7) 



This mass function explicitly breaks the Poincare symmetry of 2k + 1 dimensional 
spacetime, but preserves the Euclidian Poincare symmetry of 2k dimensional space- 
time. The fermion is also assumed to interact with 2fc-dimensional background gauge 
fields A^x^) which are independent of s. The Dirac equation may be written as: 



[0 + Td s + m(s)] *(x M , s) = 



(3.8) 



where Tj) is the lower dimension (d = 2k) covariant Dirac operator. The spinor ^ can 
be factorized as the product of a functions of s times spinors ipiXfj), 



V( Xfl , s) = }2 K(s)P+ + f n (a)P-] M**) . 

n 

satisfying the equations 

[d a + m(s)]b n (s) = n„f n (s) , 
[-d s + m(s)]f n (s) = n n b n (s) , 



p+ = 



i±r 



(3.9) 



and 



+ fi n )q n (x) = 



(3.10) 



(3.11) 



One might expect all the eigenvalues in eqn. (3.101 to satisfy \/j, n \ > 0{m), since that 
is the only scale in the problem. However, there is also a solution to eqn. (3.101 with 
eigenvalue /1 — given by 

6 = iVe~ /o" m( - s "> ds ' = ^e" m|s| . (3.12) 
This solution is localized near the defect at s — 0, falling off exponentially fast away 



from it. There is no analogous solution to eqn. (3.101 of the form 



f ~ e+fo m(s')ds' > 

since that would be exponentially growing in |s| and not normalizable. Therefore as 
seen from eqn. (3.11 1 the spectrum consists of an infinite tower of fermions satisfying 
the d — 2k Dirac equation: massive Dirac fermions with mass 0(m) and higher, plus 
a single massless right-handed chiral fermion. The massless fermion is localized at the 
defect at s = 0, whose profile in the transverse extra dimension is given by eqn. (3.12); 



the massive fermions are not localized. Because of the gap in the spectrum, at low 
energy the accessible part of the spectrum consists only of the massless RH chiral 
fermion. 
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Exercise 3.1 Construct a d = 2k + 1 theory whose low energy spectrum possesses a single 
light d = 2k Dirac fermion with mass arbitrarily lighter than the domain wall scale m. There 
is more than one way to do this. 



Some comments are in order: 

» It is not a problem that the low energy theory of a single right-handed chiral 
fermion violates parity in d = 2k since the mass for 'J breaks parity in d = 2k + 1; 

• Furthermore, nothing is special about right-handed fermions, and a left handed 
mode would have resulted if we had chosen the opposite sign for the mass in 



eqn. (3.7). This makes sense because choosing the opposite sign for the mass can 
be attained by flipping the sign of all the space coordinates: a rotation in the 
(2k + 1) dimensional theory, but a parity transformation from the point of view 
of the 2/c-dimensional fermion zeromode. 

• The fact that a chiral mode appeared at all is a consequence of the normalizability 
of exp(— J Q m(s')ds'), which in turn follows from the two limits m(±oo) being 
nonzero with opposite signs. Any function m(s) with that boundary condition 
will support a single chiral mode, although in general there may also be a number 
of very light fermions localized in regions wherever |m(s)| is small — possibly 
extremely light if m(s) crosses zero a number of times, so that there are widely 
separated defects and anti-defects. 

• Gauge boson loops will generate contributions to the fermion mass function which 
are even in s. If the coupling is sufficiently weak, it cannot effect the masslessness 
of the chiral mode. However if the gauge coupling is strong, or if the mass m 
is much below the cutoff of the theory, the radiative corrections could cause the 
fermion mass function to never change sign, and the chiral mode would not exist. 
Or it could still change sign, but become small in magnitude in places, causing 
the chiral mode to significantly delocalize. An effect like this can cause trouble 
with lattice simulations at finite volume and lattice spacing; more later. 

3.3 Domain wall fermions and the Callan-Harvey mechanism 



Now turn on the gauge fields and see how the anomaly works, following (Callan and 



Harvey, 1985). To do this, I integrate out the heavy modes in the presence of a back- 
ground gauge field. Although I will be interested in having purely 2fc-dimcnsional gauge 
fields in the theory, I will for now let them be arbitrary 2k + 1 dimensional fields. And 
since it is hard to integrate out the heavy modes exactly, I will assume perform the 
calculation as if their mass was constant, and then substitute m(s); this is not valid 
where m(s) is changing rapidly (near the domain wall) but should be adequate farther 



away. Also — in departure from the work of (Callan and Harvey, 19851, I will include 
a Pauli-Villars field with constant mass M < 0, independent of s; this is necessary to 
regulate fermion loops in the wave function renormalization for the gauge fields, for 
example. 

When one integrates out the heavy fields, one generates a Chern Simons operator 
in the effective Lagrangian, as discussed in §2.4| 
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c ™={W)\ + w\) 0cs={e{s) - 1)0cs (3 - 13) 

Note that with M < 0, the coefficient of the operator equals —2 on the side where 
m(s) is negative, and equals zero on the side where it is positive. For a background 
U(l) gauge field one finds in Euclidian spacetime: 

d = 3 : O CS = ~^ abc (A a d a A c ) , (3.14) 
e 3 

d = 5 : O CS = -^2ea6cde(A a 9 6 A c a d A e ) . (3.15) 
Differentiating Cqs by A M and dividing by e gives the particle number current: 

! — ^£abc(Fb c ) d = 3 

(3.16) 
_ 63772 £abcde(Fb c Fde) d = 5 

where I use Latin letters to denote the coordinates in 2fc + 1 dimensions, while Greek 
letters will refer to indices on the 2fc-dimcnsional defect. So when we turn on back- 
ground 2k dimensional gauge fields, particle current flows either onto or off of the 
domain wall along the transverse s direction on the left side (where m(s) = — m). If 
we had regulated with a positive mass Pauli Villars field, the current would flow on the 
right side. But in cither case, this bizarre current exactly accounts for the anomaly. 
Consider the case of a 2-dimensional domain wall embedded in 3-dimensions. If we 
turn on an E field we know that from the point of view of a 2d creature, RH Weyl 



particles are created, where from eqn. (2.14) 



We see from eqn. (3.16) this current is exactly compensated for by the Chern Simons 



(CS) 

current J2 = -f-^^i/F^ which flows onto the domain wall from the — s = —x-i side. 
The total particle current is divergenceless. 

This is encouraging: (i) we managed to obtain a fermion whose mass is zero due 
to topology and not fine tuning; (ii) the low energy theory therefore has a chiral 
symmetry even though the full 3d theory does not; (iii) the only remnant of the 
explicit chiral symmetry breaking of the full theory is the anomalous divergence of the 
chiral symmetry in the presence of gauge fields. One drawback though is the infinite 
dimension in the s direction, since we will eventually want to simulate this on a finite 
lattice; besides, it is always disturbing to see currents streaming in from s — —00! 
One solution is to work in finite (2k + 1) dimensions, in which case we end up with 
a massless RH mode stuck to the boundary on one side and a LH mode on the other 
(which is great for a vector like theory of massless Dirac fermions, but not for chiral 
gauge theories). This is what one does when simulating domain wall fermions. The 
other solution is more devious, leads to the "overlap operator" , and is the subject of 
another day's lecture. 
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3.3.1 Domain wall fermions on a slab 

To get a better understanding for how the theory works, it is useful to consider a com- 
pact extra dimension. In particular, consider the case of periodic boundary conditions 
^(x^,, s + 2so) = ^(x^s); we define the theory on the interval — sq < s < sq with 
^(x^, —so) = ^(x^, so) and mass m(s) = mA. Note the the mass function m(s) now 
has a domain wall kink at s — and an anti-kink at s = ±sq. There are now two exact 
zeromode solutions to the Dirac equation, 

b ( S ) = Ne~ I- Bo m{s')d S ' /q(s) = Ne +r^ Q rn(s')ds' 

Both solutions are normalizable since the transverse direction is finite; &o corresponds 
to a right-handed chiral fcrmion located at s = 0, and /o corresponds to a left-handed 
chiral fcrmion located at s — ±so- However, in this case the existence of exactly 
massless modes is a result of the fact that f^*° m(s) ds — which is not a topological 
condition and not robust. For example, turning on weakly coupled gauge interactions 
will cause a shift the mass by 5m(s) cx am (assuming m is the cutoff) which ruins this 
property. However: remember that to get a mass in the 2/c-dimensional defect theory, 
the RH and LH chiral modes have to couple to each other. The induced mass will be 

Sfio ~ 5m ( dsb Q (s) /o(s) = 5m N 2 ~ am x — — - = m rcs (3.19) 
J cosh|msoj 

which vanishes exponential fast as (Mso) — > oo. Nevertheless, at finite so there will 
always be some chiral symmetry breaking, in the form of a residual mass, called m res . 
If, however, one wants to work on a finite line segment in the extra dimension instead 
of a circle, we can take an asymmetric mass function, 

m(s) = \- m - S ;^° (3.20) 

I +00 < S < S 

This has the effect of excluding half the space, so that the extra dimension has bound- 
aries at s = —so and s = 0. Now even without extra interactions, one finds 

m YCS - 2me- 2mso (3.21) 

Any matrix element of a chiral symmetry violating operator will be proportional 
to the overlap of the LH and RH zeromode wave functions, which is proportional to 
m rcs . On the lattice the story of m res is more complicated — as discussed in §3.4 — 
both because of the discretization of the fermion action, and because of the presence of 
rough gauge fields. Lattice computations with domain wall fermions need to balance 
the cost of simulating a large extra dimension versus the need to make m res small 
enough to attain chiral symmetry. 

3.3.2 The (almost) chiral propagator 

Before moving to the lattice, I want to mention an illuminating calculation by Luscher 



(Luscher, 2000 a I who considered noninteracting domain wall fermions with a semi- 



infinite fifth dimension, negative fcrmion mass, and and LH Weyl fermion zeromode 
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bound to the boundary at s = 0. He computed the Green function for propagation of 
the zeromode from (x, s = 0) to (y, s = 0) and examined the chiral properties of this 
propagator. The differential operator to invert should be familiar now: 

D 5 = 04 + 75<9 S ttl , s > . (3.22) 

We wish to look at the Green function G which satisfies 

D 5 G(x,s;y,t) = S^(x-y)S(s-t) , P+G(x, 0; y, t) = . (3.23) 

The solution Liischer found for propagation along the boundary was 

G(x,s;y,t) = 2P_D~ X P + , (3.24) 

s=t=0 

where D is the peculiar looking operator 

D=[l + l5 e(H)] , H = l5 (0 A -m) = tf , e(0) = -J^= . (3.25) 

This looks pretty bizarre! Since if is hermitian, in a basis where H is diagonal, e(i?) = 
±1! But don't conclude that in this basis the operator is simply D = (1 ± 75) — you 
must remember, that in the basis where H is diagonal, e(H)js is not (by which I 
mean (m\e(H)j^\n) is in general nonzero for m ^ n in the H cigenstate basis). In fact, 



eqn. (3.25 1 looks very much like the overlap operator discovered some years earlier and 
which we will be discussing soon. 

A normal Weyl fermion in four dimensions would have a propagator P_(^ 4 ) _1 P + ; 
here we see that the domain wall fermion propagator looks like the analogous object 
arising from the fermion action , with D playing the role of the four-dimensional 
Dirac operator 4 . So what are the properties of Dl 

• For long wavelength modes (e.g. k m) we can expand D in powers of $ 4 and 
find 

D=-(0 4 -^ + ..) , (3.26) 



m \ 2m 

which is reassuring: we knew that at long wavelengths we had a garden variety 
Weyl fermion living on the boundary of the extra dimension (the factor of 1 /m is 
an unimportant normalization). 
• A massless Dirac action is chirally invariant because {75, 4 } — 0. However, the 
operator D does not satisfy this relationship, but rather: 

{75,^1 = ^75^, (3.27) 

or equivalently, 

{75,^} = 75- (3.28) 



This is the famous Ginsparg- Wilson equation, first introduced in context of the 
lattice (but not solved) many years ago (Ginsparg and Wilson, 19821. Note the 
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right hand side of the above equations encodes the violation of chiral symmetry 



that our Weyl fermion experience; the fact that the right side of eqn. (3.281 is 
local in spacetime implies that violations of chiral symmetry will be seen in Green 
functions only when operators are sitting at the same spacetime point. We know 
from our previous discussion, the only chiral symmetry violation that survives to 
low energy in the domain wall model is the anomaly, and so it must be that the 



chiral symmetry violation in eqns. ( 3.27 )-( 3.28 1 encode the anomaly and nothing 
else, at low energy]]] 

3.4 Domain wall fermions on the lattice 

The next step is to transcribe this theory onto the lattice. If you replace continuum 
derivatives with the usual lattice operator D — > g(V* + V) (where V and V* are 
the forward and backward lattice difference operators respectively) then one discov- 
ers... doublers! Not only are the chiral modes doubled in the 2k dimensions along the 
domain wall, but there are two solutions for the transverse wave function of the zero 
mode, &o(s), one of which alternates sign with every step in the s direction and which 
is a LH mode. So this ends up giving us a theory of naive fermions on the lattice, only 
in a much more complicated and expensive way! 

However, when we add Wilson terms |V*V for each of the dimensions, things 
get interesting. You can think of these as mass terms which are independent of s but 
which are dependent on the wave number k of the mode, vanishing for long wavelength. 
What happens of we add a fc-dependent spatially constant mass Am(k) to the step 



function mass m(s) — me(s)? The solution for bo(s) in eqn. (3.121 for an infinite extra 
dimension becomes 

b = Ne~ /•'M , ')+M#' , (3.29) 

which is a normalizablc zeromode solution — albeit, distorted in shape — so long as 
|Am(fe)| < m. However, for |Am(fc)| > m, the chiral mode vanishes. What happens 
to it? It becomes more and more extended in the extra dimension until it ceases to 
be normalizablc. What is going on is easier to grasp for a finite extra dimension: as 
|Am(fe)| increases with increasing k, eventually the bo zeromode solution extends to 
the opposing boundary of the extra dimension, when |Am(fc)| ~ (m — 1/so). At that 
point it can pair up with the LH mode and become heavy. 

So the idea is: add a Wilson term, with strength such that the doublers at the 
corners of the Brillouin zone have |Am(fc)| too large to support a zeromode solu- 
tion. Under separation of variables, one looks for zeromode solutions with ^(x, s) — 
e lpx tp±(s)tp± with Tip± — ±ip. One then finds (for r = 1) 

M± = 0, -4>±{sTl) + (m ee (s) + l)^ ± {s) =0 , (3.30) 

1 A lattice soluti on to eqn. J3.27t (the only solution in existence) is the overlap operator discov- 
ered by Neuberger jNeuberger ,T558a, Neuberg er, 19986>; it was a key reformulation of earlier work 
( Narayanan and Neuberger, 1993 Narayanan and Neuberger, 1995 ) on how to represent domain wall 
fermions with an infinite extra dimension (and therefore exact chiral symmetry) in terms of entirely 
lower dimensional variables. We will discuss overlap fermions and the Ginsparg- Wilson equation fur- 
ther in the next lecture. 
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Fig. 3.1 Domain wall fermions in d — 2 on the lattice: dispersion relation plotted in the 
Brillouin zone. Chiral modes exist in white regions only. For < \m/r\ < 2 there exists a 
single RH mode centered at (fa, fa) = (0,0). for 2 < \m/r\ < 4 there exist two LH modes 
centered at (fa, fa) = (7r, 0) and (fa, fa) = (0, n); for 4 < \m/r\ < 6 there exists a single RH 
mode centered at (fa, fa) = (77, 71"). For jm/r| > 6 there are no chiral mode solutions. 



where 



m e ff(s) = rae(s) + VVl - cosp^) = me(s) + F(p) 



Solutions of the form </>±(s) = z± are found with 

z± = (1 + ™ c ff(s)) Tl = (1 + me(s) + F(p))^ 



(3.31) 



(3.32) 



they are normalizable if |z| £ ( s ) < 1. Solutions are found for ip + only, and then provided 
that m is in the range F(p) < m < F(p)+2. (For r ^ 1, this region is found by replacing 
m — * m/r.) However, even though the solution is only found for ip + , the chirality of 
the solutions will alternate with corners of the Brillouin zone, just as we found for 



naive fermions, eqn. (3.5). The picture for the spectrum in 2d is shown in Fig. 3.1 



It was first shown in (Kaplan, 1992) that doublers could be eliminated for domain 



wall fermions on the lattice; the rich spectrum in Fig. 3.1 was worked out in (Jansen 
and Schmaltz, 1992[ ), where for 4d they found the number of zeromode solutions to be 
the Pascal numbers (1,4,6,4, 1) with alternating chirality, the critical values for \m/r\ 
being 0, 2, ... , 10. One implication of their work is that the Chern Simons currents 
must also change discontinuously on the lattice at these critical values of |m|/r; indeed 
that is the case, and the lattice version of the Callan-Harvey mechanism was verified 



analytically in (Golterman, Jansen and Kaplan, 1993) 



Fig. |3 - 1 1 suggests that chiral fermions will exist in two spacetime dimensions so 
long as < |m/r| < 6, with critical points at |m/r| = 0, 2, . . . , 6 where the numbers 
of massless flavors and their chiralities change discontinuously. In four spacetime di- 
mensions a similar calculation leads to chiral fermions for < |m/r| < 10 with critical 
points at \m/r\ = 0,2,..., 10. However, this reasoning ignores the gauge fields. In 
perturbation theory one would expect the bulk fermions to obtain a radiative mass 
correction of size 5m ~ O(a) in lattice units, independent of the extra dimension s. 
Extrapolating shamelessly to strong coupling, one then expects the domain wall form 
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Fig. 3.2 A sketch of the possible phase structure of QCD with Wilson fermions where the 
shaded region is the Aoki phase — pictured extending to the continuum limit (left) or not 
(right). When using Wilson fermions one attempts to tune the fermion mass to the phase 
boundary (arrow) to obtain massless pions; this is only possible in the continuum limit of the 
picture on the left is correct. For domain wall fermions chiral symmetry results at infinite L a 
when one simulates in any of the regions marked with an "X". There are six "fingers" in this 
picture instead of five due to the discretization of the fifth dimension. 

of the mass to be ruined when a ~ 1 for \m/r\ ~ (2n +1), n = 0, . . . ,4 causing a loss 
of chiral symmetry; near the critical points in |m/r| the critical gauge coupling which 
destroys chiral symmetry will be smaller. 

While qualitatively correct, this argument ignores the discrete nature of the lattice. 
On the lattice, the exponential suppression m ros ~ exp(— 2mso) found in eqn. (3.191 
is replaced by T Ls = cxp(—L s h), where T is a transfer matrix in the fifth dimension 
which is represented by L s lattice sites. Good chiral symmetry is attained when h 
exhibits a "mass gap" , i.e. when all its eigenvalues are positive and bounded away 
from zero. However one finds that at strong coupling, rough gauge fields can appear 
which give rise to near zero-modes of h, destroying chiral symmetry, with m rcs oc 1 / L s 
(Christ, 2006 Antonio et al, 20081. To avoid this problem, one needs to work at 



weaker coupling and with an improved gauge action which suppresses the appearance 
of rough gauge fields. 

At finite lattice spacing the phase diagram is expected to look something like in 
Fig. 3.2 where I have plotted m versus g 2 , the strong coupling constant. On this 
diagram, g 2 — > is the continuum limit. Domain wall fermions do not require fine 
tuning so long as the mass is in one of the regions marked by an "X" , which yield 
{1,4,6,4, 1} chiral flavors from left to right. The shaded region is a phase called the 



Aoki phase (Aoki, 19841; it is presently unclear whether the phase extends to the 



continuum limit (left side of Fig. 3.2 1 or not (right side) (Golterman, Sharpe and 
Singl eton, 2005[ ). In either case, the black arrow indicates how for Wilson fermions on 
tunes the mass from the right to the boundary of the Aoki phase to obtain massless 
pions and chiral symmetry; if the Aoki phase extends down to g 2 — than the Wilson 
program will work in the continuum limit, but not if the RH side of Fig. |3.2|pertains. 



See (Golterman and Shamir, 2000 Golterman and Shamir, 2003) for a sophisticated 



discussion of the physics behind this diagram. 

Of course, in the real world we do not see exact chiral symmetry, since quarks 
and leptons do have mass. A mass for the domain wall fermion can be included as a 
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coupling between the LH mode at s = 1, and the RH mode at s = N s 
m q [^(x, l)P+V(x, iV s ) + ^(x, N S )P_ V(x, 1)] 



(3.33) 



and correlation functions are measured by sewing together propagators from one 
boundary to itself for chiral symmetry preserving operators, or from one boundary 
to the other for operators involving a chiral flip. The latter will require insertions of 
the mass operator above to be nonzero (assuming a neglible m rcs ) — just like it should 
be in the continuum. 

3.4.1 Shamir's formulation 

Domain wall fermions are used by a number of lattice collaborations these days, using 



the formulation of Shamir (Shamir, 1993 Furman and Shamir, 1995), which is equiva- 



lent to the continuum version of domain wall fermions on a slab described above. The 
lattice action is given by: 



5 N„ 



miptp 



6=1 



s=l 



f 

-ipd b * d b ip 



(3.34) 



where the lattice coordinate on the 5d lattice is n = {x, s}, x and s being the 4d and 
fifth dimension lattice coordinates respectively. The difference operators are 



d b i>(n) = ip(n + fi b ) - -0(n) 



3 b *V>(n) = ip(n) - ip(n - fi b ) 



(3.35) 



where fib is a unit vector in the x b direction. In practice of course, these derivatives are 
gauged in the usual way by inserting gauge link variables. The boundary conditions 
are defined by setting fields to zero on sites with s = and s = N s + 1. 1 have reversed 
the sign of m and r from Shamir's original paper, since the above sign for r appears to 
be relatively standard now. For domain wall fermions, m has the opposite sign from 
standard Wilson fermions, which is physics, not convention. The above action gives 
rise to a RH chiral mode bound to the s = 1 boundary of the lattice, and a LH chiral 
mode bound at the s — N s boundary. 

3.4.2 The utility of domain wall fermions 

Theoretically, chiral symmetry can be as good a symmetry as one desires if one is close 
enough to the continuum limit (to avoid derealization of the zeromode due to large 
gauge field fluctuations) and large extra dimension. In practical simulations, the ques- 
tion is whether the residual mass term can be small enough to warrant the simulation 



cost. This was reviewed dispassionately and at length in (Sharpe, 2007), and I refer 



you to that article if you are interested in finding out the details. Currenty, domain 
wall fermions are being extensively applied to QCD; for some diverse examples from 



the past year see (Yamazaki, 2009 Chiu, Hsieh and Tseng, 2009 & Gavai and Sharma, 



2009 Torok et al., 2009 Ohta, RBC and Collaborations, 2009 Cheng et al, 2009 



Chiu et al., 2009a), and a recent overview (Jansen, 2008). Another recent application 



has been to N = 1 supersymmetric Yang-Mills theory (Giedt, Brower, Catterall, Flem- 



ing and Vranas. 2008: Giedt. Brower. Catterall. Fleming and Vranas. 2009; Endrcs. 
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2008 Endres, 2009a Endres, 2009 6) based on a domain wall formulation for Majorana 



and Vranas, 2001 ) 



fermions (Kaplan and Schmaltz, 20001 and earlier numerical work (Fleming, Kogut 



4 



Overlap fermions and the 
Ginsparg-Wilson equation 



4.1 Overlap fermions 

We have seen that the low energy limit of a domain wall fcrmion in the limit of large 

extra dimension is a single massless Dirac fermion, enjoying the full extent of the chiral 

symmetry belonging to massless fermions in the continuum. In this low energy limit, 

the effective theory is four-dimensional if the original domain wall fermion lived in five 

dimensions. One might wonder whether one could dispense with the whole machinery 

of the extra dimension and simply write down the low energy four-dimensional theory 

to start with. Furthermore, one would like a four dimensional formulation with exact 

chiral symmetry, which could only occur for domain wall fermions with infinite extent 

in the time direction, which is not very practical numerically! 

Neuberger and Narayanan found an extremely clever way to do this, leading to 

the four-dimensional "overlap operator" which describes lattice fermions with perfect 

chiral symmetry. The starting point is to consider a five dimensional fcrmion in the 

continuum with a single domain wall, and to consider the fifth dimension to be time 

(after all, it makes no difference in Euclidian space). Then 75(^4 + 777,(5)) looks like the 

Hamiltonian, where s is the new time coordinate, and m(— oo) = —mi, m(oo) = +ni2, 

where mi ; 2 > 0. The path integral projects onto ground states, and so the partition 

function for this system is Z = (CI, — mi|f2, +1712), where the state |0,m) is the ground 

state of H4,(m) = 75(4^4 + m) We know that this should describe a massless Weyl 

fcrmion. Note that the partition function is in general complex with an ill-defined 

phase (we can redefine the phase of —mi) and |f2,m,2) separately and arbitrarily). 

If we now instead imagine that the fermion mass function m(s) exhibits a wall-antiwall 

pair, with the two defects separated infinitely far apart, we recognize a system that will 

1 1 2 

have a massless Dirac fcrmion in the spectrum, and Z — \(Q, — mi |f2, +7772) , which 
is real, positive, and independent on how we chose the phase for the groundstates. 

We can immediately transcribe this to the lattice, where we replace with the 
four dimensional Wilson operator, 

H(m) = -y 5 (D w + m) = 75 - T -D\ + m) (4.1) 

with Dfj, being the symmetric covariant derivative on the lattice, and D" 1 being the 
covariant lattice Laplacian. Note that 7i(m) is Hermitian, and so its eigenvalues are 
real. Furthermore, one can show that it has equal numbers of positive and negative 
eigenvalues. 
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We can account for the 7-matrix structure of Ti{m) explicitly in a chiral basis 
where 75 = 03 ® 1: 

B + m C 



where B = — |V 2 is the Wilson operator and C — D^cr^ where <r M = {i,a}. For 
simplicity for (f2,m2| one can take m 2 — > 00, in which case TL ~ +m.27s- 

We know that Z = (ft, — mi|f2, +TO2) will represent a massless Dirac fermion on 
the lattice, so long as < mj < 2r, with 7712 arbitrary. The groundstates of interest 
may be written as Slater determinants of all the one-particle wave functions with 
negative energy. Let us designate the one-particle energy eigenstates of TL{— mi) and 
JKjn-i) to be \n,—mi) and |n, m^) respectively, with 

(n,m 2 \n',-m 1 ) = U nnl = rf\ , U*U = 1, (4.3) 

where the block structure of U is in the same 7-matrix space that we introduced 
in writing TL in block form, eqn. (4.2). Now, we want to only fill negative energy 
eigenstates, so it is convenient to introduce the sign function 

c{x)s vk- (44) 



With mi — > 00 we have 



- - \\\ \ ^ 



Assuming 7i(— mi) has no exact zeromodes then, it follows that all eigenvalues come 
in ± pairs (just like the operator 75) and we can choose our basis \n, —mi) so that 

e(H(m 2 )) = C/ 75 C/t =U Q °^ [/t . (4 . 6) 

Therefore the Slater determinant we want is 

Z = I (f2, jri2|f2, — mi) | 2 

= j det C/ 22 1 2 
= det ,5 t det 6 

= det fl + ^y-mOn . (4.7) 



Some steps have been omitted from ths derivation (Narayanan, 2001 1; see exercise 4.2 



Exercise 4.1 Prove the assertion that if Tri{— mi) has no zeromodes, it has equal numbers 
of positive and negative eigenvalues. 
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Exercise 4.2 You should prove the last step in eqn. (4.7|, breaking it down to the following 
steps: 

(a) Show that det 5 ' = det a det ; 

(b) ...so that det<5 f det S = det S det a det LT* = det [~{U + -y 5 [J'y 5 )U'<] ; 



(c) ...which combines with eqn. ( 4.6 1 to yield eqn. ( 4.7 1 



On the other hand, Z oc det D, where D is the fermion operator. So we arrive at 
the overlap operator (dropping the subscript from mi): 

D = l + 75 £(W(-m)) 

= 1 + 7 H( ~ m) 

= 1+ - ^ w ~ — . (4.8) 

V (Ao - m^CAu - m) 

a remarkable result. It was subsequently shown explicitly that this fermion operator 
can be derived directly from lattice domain wall fermions at infinite wall separation 



(Neuberger, 1998c Kikukawa and Noguchi, 1999). Recall from our discussion of do- 



main wall fermions that at least for weak gauge fields, we need < m < 2r in order 
to obtain one flavor of massless Dirac fermion (where I have set the lattice spacing 
a = 1). 

4.1.1 Eigenvalues of the overlap operator 

Recall that the eigenvalues of the Dirac operator in the continuum are ±?A„ for real 
nonzero A„, plus n + RH and n_ LH zero modes, where the difference is constrained 
by the index theorem to equal the topological winding number of the gauge field. Thus 
the spectrum looks like a line on the imaginary axis. What does the spectrum of the 
overlap operator look like? Consider 

(D-l) t (D-l) = e(W) 2 = l . (4.9) 

Thus (D — 1) is a unitary matrix and the eigenvalues of D are constrained to lie on a 
circle of unit radius in the complex plane, with the center of the circle at z = 1. If you 
put the lattice spacing back into the problem, D — > aD in the above expression to get 
the dimensions right, and so the eigenvalues sit on a circle of radius 1/a centered at 
1/a. Thus, as a — > the circle gets bigger, and the eigenvalues with small magnitude 
almost lie on the imaginary axis, like the continuum eigenvalues. See the problem 
below, where you are to show that the cigcnfunctions of D with real eigenvalue arc 
chiral. 

4.1.2 Locality of the overlap operator 



If just presented with the overlap operator eqn. (4.7 1 without knowing how it was 



derived, one might worry that its unusual structure could entail momentum space 
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singularities corresponding to unacceptable nonlocal behavior in coordinate space. 
(From its derivation from domain wall fermions this would be very surprising for 
sufficiently weakly coupled gauge fields, since the domain wall theory looks well defined 
and local with a mass gap.) The locality of the overlap operator (i.e. that it falls 



off exponentially in coordinate space) was proven analytically in (Hernandez, Jansen 
and Luscher, 1999) , under the assumption of sufficiently smooth gauge link variables, 



namely that |1 — U\ < 1/30. They also claimed numerical evidence for locality that 
was less restrictive. 

4.1.3 The value of m and the number of fermions 

For domain wall fermions we found the interesting phase structure as a function of m/r, 
where in the intervals between the critical values m/r — {0, —2, —4 . . . , — 2(2fc + 1)} 
there were {1, 2k, . . . , 1} copies of chiral fermions with alternating chirality, as shown 
in Fig. |3.1| One would expect something analogous then for the overlap operator, 
since it is equivalent to a domain wall fermion on a 2/c-dimcnsional lattice with infinite 
continuous extra dimension. The equation of motion for the domain wall modes is 
slightly different than found in eqn. ( |3.30[ ) due to the continuous dimension: 



jW± = 0, ±(t>' ± (s) + (m eS (s) + l)ct>±(s) = , (4.10) 



where 



m e s(s) = me(s) + r V^(l — cosp M ) = me(s) + rF(p) . (4-11) 
i" 

Solutions are of the form (f)±(s) — e T f ro e f f (t)dt_ with r F(p) > 0, 0_ is never nor- 
malizable, while for <j) + , normalizability requires m/r > F(p). Thus the overlap op- 
erator eqn. ( |4.8[ ) represents {1, 1 + 2k, . . . ,2 2k } massless Dirac fermions for m/r in 
the intervals (0, 2), (2, 4), ... , (4k, oo). In 2k = 4 dimensions, these flavor numbers are 
{1,5,11,15,16}. 

4.1.4 Simulating the overlap operator 

The overlap operator has exact chiral symmetry, in the sense that it is an exact solution 
to the Ginsparg Wilson relation, which cannot be said for domain wall fermions at finite 
iV s ; furthermore, it is a four-dimensional operator, which would seem to be easier to 
simulate than a 5d theory. However, the inverse square root of an operator is expensive 
to compute, and requires some approximations. The algorithms for computing it are 
described in detail in an excellent review by A. Kennedy (Kennedy, 2006). Amusingly, 



he explains that the method for computing the overlap operator can be viewed as 
simulating a five-dimensional theory, albeit one with more general structure than the 
domain wall theory. For a recent review comparing the computational costs of different 



lattice fermions, see the recent review (Jansen, 20081. 



Exercise 4.3 Show that the overlap operator in eqn. (4.81 has the following properties: 

(a) At zero gauge field and acting on long wavelength fermion modes, D ~ $ A , the ordinary 
Dirac operator for a massless fermion. 
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(b) It satisfies the Ginsparg-Wilson equation, eqn. (3.271: 

{j 5 ,D} = D l5 D . 



(4.12) 



Exercise 4.4 

(a) Show that one can write D = 1 + V where V^V = 1, and that therefore D can be diagonalized 
by a unitary transformation, with its eigenvalues lying on the circle 2 = 1 + e . 

(b) Show that, despite D being non-hermitian, normalized eigenstates satisfying D\z) = z\z) with 
different eigenvalues are orthogonal, satisfying (z'\z) — 8 z i z 

(c) Show that if D\z) = z\z) then D^\z) = z*\z) 

(d) Assuming that 75D75 = D\ show that (z^slz) = unless z = or z = 2, in which case 
(z\ l5 \z) = ±1 



4.2 The Ginsparg-Wilson equation and its consequences 

In 1982 Paul Ginsparg and Kenneth Wilson wrote a paper about chiral lattice fermions 
which was immediately almost completely forgotten, accruing 10 citations in the first 
ten years and none in the subsequent five; today it is marching toward 700 citations. 
The reason for this peculiar history is that they wrote down an equation they spec- 
ulated should be obeyed by a fermion operator in the fixed point action of a theory 
tuned to the chiral point — but they did not solve it. After domain wall and over- 
lap fermions were discovered in the early 1990s, it was realized that they provided a 
solution to this equation (the domain wall solution only being exact in the limit of 
infinite extra dimension). Shortly afterward, M. Liischer elaborated on how the salient 
features of chirality flowed from the Ginsparg-Wilson equation — in particular, how 
anomalies and multiplicative mass renormalization were consequences of the equation, 
which provided a completely explicit four-dimensional explanation for the success of 
the overlap and domain wall fermions. 

4.2.1 Motivation 

A free Wilson fermion with its mass tuned to the critical value describes a chiral 
fermion in the continuum. As we have seen, chiral symmetry does not exist on the 
lattice, but its violation is not evident at low energy, except through correctly repro- 
ducing the anomaly. However, imagine studying this low energy effective theory by 
repeatedly performing block spin averages. One would eventually have a lattice theory 
with all the properties one would desire: chiral fermions and chiral anomalies. What is 
the fermion operator in this low energy theory, and how does it realize chiral symme- 
try? Motivated by this question, Ginsparg and Wilson performed a somewhat simpler 
calculation: they took a continuum theory with chiral symmetry and anomalies, and 
performed a average of spacetime cells to create a lattice theory, and asked how the 
chiral symmetry in the original theory was expressed in the resulting lattice theory. 
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The starting point is the continuum theory 

Z= f [dip] [#]e" s ^ (4.13) 



I assume there are Nf identical flavors of fermions, and that S is invariant under the 
full U(Nf) x U(Nf) chiral symmetry. We define tp n to be localized averages of tp, 

tp n = J d 4 xiP(x)f(x- an) (4.14) 

where /(x) is some function with support in the region of |x| < a. Then up to an 
irrelevant normalization, we can rewrite 

J J n 

= [l[dXn dxn e- s ^"^ = e-* D * , (4.15) 

J n 

where a is a dimensionful parameter, where D is the resulting lattice fermion operator. 
Since there are Nf copies of all the fields, the operator D is invariant under the vector 
U(Nf) symmetry, so that if T is a U(Nf) generator, [T, D] = 0. The lattice action is 
therefore defined as 

* Dx = { [dip] [dtp] e~[ E " «(x»-?-)(x»-iW-S(iM)] ; (4.i 6 ) 



e 

Note that explicit chiral symmetry breaking has crept into our definition of Si a t 
through the fermion bilinear we have introduced in the Gaussian in order to change 
variables. 

Now consider a chiral transformation on the lattice variables, Xn - ► e l£75T x n , 
Xn — * Xne" 75T , where T is a generator for a U(Nf) flavor transformation. This is 
accompanied by a corresponding change of integration variables ip, ip: 

e -Xe leJ5T De^ T X = J ^ ^ j / eA Tr T g - [E„ a (x n -^ n )e 2 " TsT (x„-0„)-.S(,/>,$)] _ 

(4.17) 

where .4 is the anomaly due to the non-invariance of the measure [dtp] [dtp] as computed 
by Fujikawa (Fujikawa, 19791: 

A = j^eep-yS Tr F aP F l5 (4.18) 

with 

' A=2u , (J.IOi 
^ being the topological charge of the gauge field. 
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Expanding to linear order in e gives 



X{ 7 5, D}T X e-* D x = J [#] [#] ^2v Tr T + £ [(Xn - ^n)2a 75 T( X „ - n )] 



x exp 
^ ( -h> Tr I' 



m 

2 (5 



0m) (X 

<5 



75 T 



-0m) - S(lp, 1p) 
-XD X 

-XD X 



Tr l5 DT + 2vTvT- -x n D l5 DT Xn e 
a 



Defining a = 2/ a this yields the operator identity 

({ 75 , D} — a D l5 D) T = (Tr l5 DT + 2v Tr T) 



(4.20) 



(4.21) 



If T is taken to be a traceless generator of U(Nf), multiplying both sides by T and 
taking the trace yields the Ginsparg-Wilson equation: 



{ 75 ,D} = aD l5 D 



(4.22) 



If on the other hand we take T to be the unit matrix and use eqn. (4.221 we find 

Tr l5 D = 2N f v , (4.23) 
This latter equation was not derived in the original Ginsparg-Wilson paper; from our 



discussion of the index theorem eqn. (2.23), it follows from eqn. (4.231 that Tr 75.D = 
(n_ — n + ), where n± are the number of ± chirality zeromodes. We will see that that 
is indeed the case. 



Note that the GW relation eqn. (4.221 is the same equation satisfied by the over- 



lap operator (Neuberger, 1998 & I — and therefore by the domain wall propagator at 



infinite wall separation on the lattice, being equivalent as shown in (Neuberger, 1998a 



Neuberger, 1998 c) — as well as by the infinitely separated domain wall propagator 



in the continuum (Luscher, 2000a). In fact, the general overlap operator derived by 
Neuberger 

D = 1 + 75 e(W) (4.24) 
is the only explicit solution to the GW equations that is known. 
4.2.2 Exact lattice chiral symmetry 

Missing from the discussion so far is how the overlap operator is able to ensure mul- 
tiplicative renormalization of fermion masses (and similarly, multiplicative renormal- 
ization of pion masses). In the continuum, both phenomena follow from the fact that 
fermion masses are the only operators breaking an otherwise good symmetry. The 
GW relation states exactly how chiral symmetry is broken on the lattice, but does 
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not specify a symmetry that is exact on the lattice and capable of protecting fermion 
masses from additive renormalization. 

Luscher was able to solve this problem by discovering the GW relation implied the 
existence of an exact symmetry of the lattice action: J tpDip is invariant under the 
transformation 

H = 75 (l + \D) ip, Si> = $ (l - |d) 75 . (4.25) 

Note that this becomes ordinary chiral symmetry in the a — > limit, and that it is 
broken explicitly by a mass term for the fermions. 

4.2.3 Anomaly 

If this symmetry were an exact symmetry of the path integral, we would run afoul of 
all the arguments we have made so far: it becomes the anomalous U(1)a symmetry in 
the continuum, so it cannot be an exact symmetry on the lattice! The answer is that 
this lattice chiral transformation is not a symmetry of the measure of the lattice path 
integral: 



= [#][#] (Tr [75 ( 



-D 



Tr 



1- 2^' 75 



= [#][#] x aTr 7 5 D , (4.26) 

where I used the relation ddet M/dx = det[Af] Tr M~ 1 dM / dx. Unlike the tricky non- 
invariance of the fermion measure in the continuum under a U(1)a transformation - 
which only appears when the measure is properly regulated — here we have a perfectly 
ordinary integration measure and a transformation that gives rise to a Jacobean with 
a nontrivial phase (unless, of course, Tr 75-D = 0). To make sense, Tr 75D must map 
into the continuum anomaly. .and we have already seen that it does, from eqn. (4.23). 



What remains is to prove the index theorem (Hasenfratz, Lalicna and Nicdcrmaycr, 



1998 Luscher, 1998), the lattice equivalent of eqn. (2.221. From exercise 4.4 it follows 



that for states \z) satisfying D \ z) — z\z) 



2N f {nf -r/l 2) ) , 



(4.27) 



(2) 

where n ± are the number of positive and negative chirality states with eigenvalue 
z — 2. We also know that 

,( 2 ) „Vh 



= Tr 75 = ^(z|7 5 |z) = (n + - n_) + (n^> - vT') , 



(4.28) 



where n± are the number of ± chirality zeromodes at z = 0. Therefore we can write 

Tr 75D = 2(n_ - n+) . (4.29) 



Substituting into eqn. (4.23) we arrive at the lattice index theorem, 

(n_ — n + ) = i>Nf 



(4.30) 



which is equivalent to the continuum result eqn. (2.23), and provides an interesting 



definition for the topological charge of a lattice gauge field. A desirable feature of 
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the overlap operator is the existence of exact zeromode solutions in the presence of 
topology; it is also a curse for realistic simulations, since the zeromodes make it difficult 
to sample different global gauge topologies. And while it cannot matter what the global 
topology of the Universe is, fixing the topology in a lattice QCD simulation gives rise 



to spurious effects which only vanish with a power of the volume (Edwards, 2002) 



4.3 Chiral gauge theories: the challenge 

Chiral fermions on the lattice make an interesting story whose, final chapter on chiral 
gauge theories has barely been begun. It is a story that is both theoretically amusing 
and of practical importance, given the big role chiral symmetry plays in the standard 
model. I have tried to stress that the understanding of anomalies has been the key 
to both understanding the puzzling doubling problem and its resolution. In terms of 
practical application they are more expensive than other fermion formulations, but 
have advantages when studying physics where chirality plays an important role. For a 



recent review comparing different fermion formulations, see ( Jansen, 2008 1 . 

While domain wall and overlap fermions provide a way to represent any global 
chiral symmetry without fine tuning, it may be possible to attain these symmetries 
by fine tuning in theories with either staggered or Wilson fermions. In contrast, there 
is currently no practical way to regulate general nonabelian chiral gauge theories on 
the lattice. (There has been a lot of papers in this area, however, in the context 
of domain wall - overlap - Ginsparg-Wilson fermions; for a necessarily incomplete 



list of references that gives you a flavor of the work in this direction, see (Kaplan, 
1992; Kaplan, 1993; Narayanan and Neuberger, 1993 Narayanan and Neubergcr, 1995; 



Narayanan and Neuberger, 1996; Kaplan and Schma 



tz, 1996| |Luscher, 1999] |Aoyama| 



and Kikukawa, 1999 Luscher, 2000 &| |Kikukawa and Nakayama, 2001||Kikukawa, 2002 



Kadoh and Kikukawa, 2008| |Hasenfratz and von Allmen, 2008 1 ). Thus we lack of 



tia 
12; 



a nonperturbative regulator for the Standard Model — but then again, we think 
perturbation theory suffices for understanding the Standard Model in the real world. 
If a solution to putting chiral gauge theories on the lattice proves to be a complicated 
and not especially enlightening enterprise, then it probably is not worth the effort 
(unless the LHC finds evidence for a strongly coupled chiral gauge theory!). However, 
if there is a compelling and physical route to such theories, that would undoubtedly 
be very interesting. 

Even if eventually a lattice formulation of the Standard Model is achieved, we must 
be ready to address the sign problem associated with the phase of the fermion deter- 
minant in such theories. A sign problem has for years plagued attempts to compute 
properties of QCD at finite baryon chemical potential; the same physics is responsible 
for poor signal/noise ratio experienced when measuring correlators in multi-baryon 
states. To date there have not been any solutions which solve this problem. We can at 
least take solace in the fact that the sign problems encountered in chiral gauge theo- 
ries and in QCD at finite baryon density are not independent! After all, the standard 
model at fixed nontrivial SU(2) topology with a large winding number can describe a 
transition from the QCD vacuum to a world full of iron atoms and neutrinos! 
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